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ANALYTICAL AND EXPERIMENTAL VIBRATION STUDIES 
OF A 1/8-SCALE SHUTTLE ORBITER 


Larry D. Pinson, Coordinator 
Langley Research Center 

SUMMARY 

Theoretical and experimental results for vibration of a l/8-scale shuttle orbiter 
model are reported. Natural frequencies and mode shapes for four symmetric vibra- 
tion modes and four antisymmetric modes are compared with predictions based on 
NASTRAN finite-element analyses. Initial predictions gave poor agreement with test 
data; an intensive investigation revealed that the major factors influencing agreement 
were out-of -plane imperfections in fuselage panels and a soft fin-fuselage connection. 
Computations with a more refined analysis indicated satisfactory frequency predictions 
for all modes studied, within 11 percent of experimental values. 

INTRODUCTION 

Knowledge of the structural dynamic characteristics of launch vehicles is funda- 
mental to the prediction of such dynamics problems as pogo, control, flutter stability, 
and the response to gusts and staging transients. Early in the development cycle, 
dynamic characteristics usually are defined by analysis, but verification by full-scale 
tests is not possible due to lack of prototype hardware. Confidence in analytical pro- 
cedures may be established early in the design through the use of subscale structural 
models, and, in addition, full-scale dynamic behavior may be approximated through 
application of principles of similitude. The concept was developed by the NASA Langley 
Research Center (LRC) using a 1/5 -scale model of the Saturn I vehicle (refs. 1 to 5) 
and was applied extensively in the development of a large Air Force launch system 
(refs. 6 to 8). Early data from a 1/10-scale replica model were used to improve ana- 
lytical methods applicable to the Apollo/Saturn V launch vehicle (refs. 9 to 16) and to 
approximate full-scale behavior. The 1/10-scale Apollo/Saturn V model subsequently 
was used as a troubleshooting tool when anomalous behavior occurred in an early flight 
(ref. 17). 

During the technology development phase, the space shuttle was recognized as 
having more complex structural dynamic characteristics than previous launch vehicles 
because four separate large elements are joined asymmetrically at a few discrete 



interfaces. This multielement configuration will have high modal density at low frequen- 
cies, and the dynamic response will have a high degree of directional coupling. Thus, 
early verification of the ability to analyze adequately the various vehicle elements and 
subsequently to couple the element characteristics for prediction of total vehicle char- 
acteristics is of great importance. In addition, if confidence in this ability is sufficient, 
a reduction in full scale testing could be effected with consequent cost reduction. Studies 
of a greatly simplified multielement configuration are reported in references 18 and 19. 

To provide a better assessment of analytical modeling procedures and data with 
which to understand the dynamic behavior of shuttlelike configurations, a l/8-scale 
dynamic model of an early shuttle configuration was built for structural dynamics inves- 
tigations at the Langley Research Center. Because of the lack of detailed definition at 
the time of construction, only overall dimensions were scaled according to principles of 
similitude. Stock sizes of materials were used to approximate varying thickness and 
nonstandard scaled dimensions. The model was designed by the Grumman Aerospace 
Corporation (GAC) under contract to LRC; subsequently, GAC built the solid rocket 
boosters and the external tank. The orbiter element was constructed at LRC. 

The purpose of the present paper is to report the results of investigations of the 
vibration behavior of the l/8-scale orbiter model. During initial vibration tests of the 
orbiter model, natural frequencies were found to be significantly different from those 
predicted by a preliminary NASTRAN finite-element model. For example, for the first 
mode the analytical natural frequency was 22 percent higher than the test natural fre- 
quency. Detailed studies of the orbiter were conducted to understand and reconcile 
these differences. The studies included free vibration tests, vibration tests under static 
preload, static -load deflection tests, and detailed analytical and experimental investiga- 
tions of structural panels and joints. Results of the investigation are presented in this 
report. 

The work reported herein was accomplished through a team effort consisting of 
several LRC, GAC, and Rockwell International Corporation (RIC) personnel. On-site 
support was provided by the Lyndon B. Johnson Space Center through RIC. Thus, vari- 
ous sections of this report have different authors and, where appropriate, the author- 
ship is so designated. These authors are Leadbetter (LRC), Sewall (LRC), Blanchard 
(LRC), Flynn (GAC), Durling (LRC), Miller (GAC), Bernstein (GAC), Housner (LRC), 
Stein (LRC), Barrett (RIC), Herr (LRC), and Robinson (LRC). In addition, off-site ana- 
lytical support was provided by Philip W. Mason, J. Zalesak, and H. G. Harris of GAC. 
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SYMBOLS AND ABBREVIATIONS 


Measxirements and calculations were made in the U.S. Customary Units. They are 
presented herein in the International System of Units (SI) with the equivalent values given 
parenthetically in the U.S. Customary Units. The followir^ symbols apply to all parts of 
the report except appendix D, which has its own list of symbols. 

A effective skin area (see eq. (E13)) 

Ab,Ac,At effective skin areas in direct stress for fuselage along bottom of 
sidewall, orbiter center line, and top of sidewall, respectively 
(see sketch (El)) 

Ap cross-sectional area of cabin-to-payload bay splice plate (see eq. (C16)) 

A]^,A 2 ,A 3 ,A 4 cross-sectional areas of fuselage longerons in sketch (El) 

Aj,A 2 cross-sectional areas of wing-spar caps in sketch (E3) 

a moment arm for simulated fin-root joint flexibility (see eqs. (C8)) 

B,C,D half- width of fuselage, keel height, and height of fuselage sidewall, 

respectively (see sketch (El)) 

b fin-root width (see sketch (C4)) 

bp length of cabin-to-payload bay lower splice plate (appendix C) 

c/c dampii^ coefficient (see table II) 

E Young's modulus 

f circular frequency, co/ 27 t, Hz 

G shear modulus 

h fin height (see sketch (C4)) 
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Ia flexural moment of inertia of fuselage or wii^ cross section about neutral 

surface based on initial analysis (see eq. (E2)) 

Ia X total area moment of inertia of wii^ cross section about neutral surface 

(parallel to xy-plane) of wing (see eq. (ElO)) 

If area moment of inertia of fin-fuselage spar clip about its cross-sectional 

neutral axis (see appendix C) 

Ix required area moment of inertia of wing cross section about neutral surface 

(parallel to xy-plane) for analytical model revised by static test data 
(see eq. (ElO) and sketch (E3)) 

ly required area moment of inertia in pitch of fuselage about neutral axis for 

analytical model revised by static test data (see eq. (E5) and sketch (El)) 

I 2 required area moment of inertm m yaw of fuselage about vertical axis 

through orbiter center line for analytical model revised by static test 
data (see eq. (E6) and sketch (El)) 

Ic area moment of inertia of wing spar caps about neutral surface (parallel 

to xy-plane) of wing cross section (see eq. (Ell)) 

Ip area moment of inertia of cabin-to -payload bay splice plate about its cross- 

sectional neutral axis (see eq. (C17)) 

Ig area moment of inertia of reduced skin for wing (see eq. (Ell)) 

[l^ stiffness matrix (see eqs. (1) and (2)) 

Kfi,Kf 2 >Kf 3 , 1 spring constants representing fin-root stiffness (see fig. 24 

Kpi f,Kfi a ) appendix C) 

Ku,Kj spring constants representing cabin-to -payload bay joint stiffness 
(see fig. 25 and eqs. (C20) and (C21)) 

^1>^2 spar -clip length and splice-plate length, respectively (see sketch (C2)) 


mass matrix (see eq. (1)) 



Mj,M 2 bendir^ moments in fin-root spar clip (see eqs. (Cl) to (C4)) 

P load 

{p} load vector (see eq. (2)) 

R reaction (see sketch (C2)) 

% j^sjtu skin thicknesses of lower cover panels, fuselage sidewall, and wing upper 
cover panels, respectively (see table XVII) and sketch (E3)) 

U strain energy defined by equation (Cl) 

(u) displacement vector (see eqs. (1) and (2)) 

W total effective skin width, E^i (see eq. (E12)) 

Wp full skin width (see sketch (E3)) 

wj effective skin width in direct stress for ith spar cap (see sketch (E3)) 

Wq displacement at center of panel (see appendix E) 

X,Y,Z Cartesian coordinates 

Xj^,X 2 axial coordinates in fin-root joint analysis (see sketch (C2)) 

z^,z^ vertical distances between wii^ neutral surface and midplanes of lower 

and upper covers of wing, respectively (see sketch (E3)) 

z vertical distance from bottom of fuselage to neutral axis (see eq. (E4)) 

a angle between spar -clip plates or splice plates (see sketch (C2)) 

^ supplement to a for splice plates (see sketch (C5)) 

^ - ®e,s " ®a,s 
6 static deflection 


5 



6^ maximum static deflection from original analysis (see eq. (E2)) 

Sq- maximum static test deflection (see eq. (E2)) 

®a,s>®e,s analytical and measured fin fore-and-aft static deflections (see sketch (C3)) 
strain due to static loading (see figs. 17 and 18) 
eigenvalue (see eq. (1)) 
eigenvector (or mode shape) 
angular frequency, 27rf, radius per second 

Subscripts: 

dyn dynamic (see fig. 24) 

m model (see table I) 

p prototype (see table I) 

Abbreviations: 

RSI rod-shear panel model with fuselage and wing carry-through skins 

33 -percent effective in direct stress, fully effective in shear; fully 
effective wing skins; joint flexibilities in fin root (model FI, fig. 24) 
and cabin-to payload bay joint (fig. 25) 

RS2 rod-shear panel model with fuselage and wing carry-through skins 

46 -percent effective in direct stress; wing skins 80-percent effective 
in direct stress; all skins fully effective in shear; same joint flexi- 
bilities as RSI 

RS3 rod-shear panel model with wii^ skin 80-percent effective in direct stress 

and shear; original fuselage (i.e., fully effective skins) and original joints 


e 

X 

(<^> 

CO 
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RS4 

RS5 

RS6 

RMl 

RM2 

F1,F2 

GAC 

LRC 

NASA 

RIC 

CPB 

CPU 

MPC 

OPS 

SPC 

TPS 

s.s. 


rod-shear panel model of same direct-stress effectiveness as RS2 with fuse- 
lage skins 46 -percent and wing skins 80 -percent effective in shear, model F2 
fin-root joint flexibility, and same cabin -to -payload bay joint flexibility 

rod -shear panel model same as RS4 with both fuselage and wing skins 
80-percent effective in shear 

rod- shear panel model same as model II of reference 21 

reduced-membrane model with all exterior fuselage and wing skins 
60 -percent effective in direct stress and shear with same joint flexi- 
bilities as in RS4 and RS5 except for model F3 fin -root joint flexibility 
replacing F2 in the antisymmetric model 

reduced -membrane model same as RMl with membrane panels in door, 
also 60-percent effective in direct stress and shear 

F3 fin-fuselage interface models shown in figure 24 
Grumman Aerospace Corporation 
Langley Research Center 

National Aeronautics and Space Administration 

Rockwell International Corporation 

cabin-to -payload bay 

central processing unit 

multiple -point constraints 

orbiting propulsion system 

single -point constraints 

thermal protection system 

simply supported (see fig. 50) 
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ORBITER TEST MODEL 


Sumner A. Leadbetter and Ulysse J. Blanchard 

This section presents results of vibration tests of the orbiter element of the 
l/8-scale space shuttle dynamic model. The tests were conducted to provide data for 
the verification of the NASTRAN finite -element representation of the subscale orbiter 
model. A description of the experimental model, description of apparatus and pro- 
cedures used to conduct tests, and a discussion of results are presented. 

Model Description 

The general arrangement of the l/8-scale shuttle model is shown in figure 1. The 
design was based on a GAC configuration developed in the space shuttle studies for 
NASA in 1972. The parallel-burn configuration had a gross lift-off weight of 21.35 MN 
(4.8 X 10® lb) and was 55.5 m (182 ft) long. The complete model consists of four major 
elements: the orbiter, external tank, and two solid rocket boosters. The orbiter model, 
which was constructed at LRC and is the subject of the present investigation, is shown 
in figure 2 without the cargo -bay door and in figure 3 with the door installed. Figure 4 
is a drawing of the general arrangement of the orbiter model, and figure 5 shows the 
fuselage under construction. The orbiter model is primarily a riveted 2024-aluminum 
structure consisting of thin nontapered skins over supporting frames and longerons and 
with minimum use of intermediate skin stiffeners. The bottom and sides of the fuse- 
lage are flat, as are surfaces of the fin stub and wing. In simplifying the design, a major 
objective was to keep the model fabrication cost within target while retaining as many of 
the significant structural dynamic characteristics as possible. Hence, only the general 
characteristics of the major orbiter components were simulated. No attempt was made 
to model local details. A more detailed description of the model is presented in refer- 
ences 20 and 21. 

Scaling . - The scaling relationships that must exist between the model and the proto- 
type are shown in table I, and they are based upon the determination of the major factors, 
presented in reference 22, which influence the response of the prototype being studied. 
These relationships follow directly from a dimensional analysis of the various parameters 
that influence the dynamic behavior of the structure and from the choice of the model 
material. Extrapolating prototype behavior from model test data is accomplished by 
directly using these scaling relationships. It should be noted, however, that because of 
design expediency some of the scaling rules have been compromised. For example, the 
local skin stiffness on the model is less than the required scaled value of the prototype 
for preventing buckling. Some liberty has also been taken in modelii^ the stiffness 
characteristics by some necessary lumpii^ in order to avoid the large expense of exact 
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scaling of very small dimensions. If local stiffeners (e.g., stringers) are completely 
eliminated, then the local stiffness of the skin is not duplicated and premature buckling 
will occtir. 

While accurate modeling of the prototype was desirable for extrapolating basic 
shuttlelike dynamic behavior, a prime object of the study was the correlation of a 
NASTRAN dynamic analysis with model test data. A complete static and dynamic 
analysis was made using NASTRAN with the structure modeled to a degree of refine- 
ment considered sufficient for preliminary design purposes. Thus, the need for direct 
scaling of the prototype design to obtain an exact model in every detail was not consid- 
ered to be crucial. Also, the shuttle design was still in a state of flux at the beginning 
of this study; hence, any attempt to model the then current vehicle exactly would not be 
greatly beneficial to the shuttle project. 

Fuselage .- The fuselage is 3.543 m (11.625 ft) long, contains 21 frame stations, 
and is constructed of 2024 aluminum. The bottom skin of the fuselage is 0.635 mm 
(0.025 in.) thick, and the sidewall and top skins are 0.508 mm (0.020 in.) thick. The 
fuselage frames in the region of the cargo bay (fig. 5) are constructed of aluminum 
sheet that has been bent to form channel sections. The tapered sidewall channel sec- 
tion and the lower portion are attached back-to-back to form a U-shaped frame. The 
major bulkheads located in the forward and aft sections of the fuselage are of stiff- 
ened sheet construction. 

Fin .- The model fin structure, which represents only the structure from the fuse- 
lage to the center of gravity of the orbiter fin, contains three spars and a closure rib at 
the tip. The webs are 0.8128 mm (0.032 in.) thick while the covers are 0.5080 mm 
(0.020 in.) thick. During tests, a 26.69-N (6-lb) ballast weight representing the remain- 
ing fin structural weight was attached to the tip closure rib (location indicated in fig. 2). 
The fin structure is riveted onto the fuselage structvire. Simple center-line clips con- 
nect the fore and aft fin spars to fuselage frames, and the center fin spar is connected 
by a clip and gusset combination at each side of the fin. 

Cargo -bay door .- The fuselage cargo-bay door (fig. 3) is comprised of seven seg- 
ments of 0.4064-mm (0.016-in.) aluminum sheet attached to semicircular frames to form 
a semicylindrical shell. The design of the eight door frames (V-shaped angles) allows 
expansion in the longitudinal direction, thereby offering little resistance to fuselage bend- 
ing but allowing the fuselage-door combination to resist torsion. When attached to the 
fuselage, the model door simulates structural properties of a prototype door in a closed 
and locked position. The door has four key stone -shaped straps riveted to each side lon- 
geron which connect to the fuselage upper longeron by means of a single screw fastener 
each, as shown in figure 6. Five fasteners attach the door skin to the fuselage at each 
end of the door. 
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Wings . - The delta wings shown under construction in figtire 7 consist of six spars 
and four ribs that are formed from 0.8128-mm (0.032-in.) 2024 aluminum sheet. The 
covers are 0.508 mm (0.020 in.) thick. The proper scaled wir^ weight, includii^ items 
such as thermal protection system (TPS) panels and landing gear, is simulated by ballast 
installations shown in the figure. 

Ballast .- Locations of other concentrated mass representations are indicated in 
figure 2. These masses include the orbiting propulsion system (OPS) ballast attachment, 
the crew -cabin area ballast, and the cargo-bay payload ballast. The payload simulation 
shown in figure 2 is representative of an intermediate-sized cargo; however, the present 
vibration investigation was conducted using a maximum payload simulation unit which 
extends the entire length of the cargo bay and represents a prototype weight of 289 kN 
(65 000 lb). The model payload is a box beam with brackets simulating orbiter payload 
attachment points. Appropriate ballast is attached to the beam to provide desired scaled 
payload mass. 


Apparatus and Test Procedure 

Suspension system . - The orbiter model was suspended in an inverted horizontal 
attitude during most of the vibration tests, as shown in figure 8. Steel cables with in-line 
soft coil springs were attached at the two fuselage interstage fittings and provided essen- 
tially free -free constraints. Rigid-body natural frequencies on this suspension system 
(f < 1.0 Hz) were much lower than the lowest structural frequency (f = 43 Hz) for all tests. 
Symmetric and antisymmetric resonances were excited by applying appropriate shaker 
forces, as shown in figure 8, for example. 

Limited vibration tests also were conducted with the model suspended vertically 
as illustrated by the sketch in figure 9. During these tests, tension loads were applied 
to the model fuselage by means of a shock cord arrangement attached to the lower sur- 
face skin and the two lower main engine fittings. This test setup was used to study the 
effect of fuselage -panel stresses on resonant frequencies. The model was excited in 
the z -direction, xz -plane, while under various tension preloads ranging from 0 to 
2670 N (600 lb). 

Instrumentation . - The instrumentation of the 1/8-scale orbiter model was designed 
to provide continuous electronic signals to define the dynamic response of the model to 
sinusoidal force inputs. The instrumentation provided a means of monitoring exciter 
force inputs and transducer (acceleration) output signals, recording these signals for 
subsequent in-depth analysis, and performing on-site data analysis. The locations of 
fixed transducers, piezoelectric accelerometers, are shown in figure 10. In general, 
the accelerometers were mounted at the lower corners of the fuselage frames (stations), 
at available stiff locations such as the juncture of wing spars and ribs, and on ballast 
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masses. At all fuselage instrument locations, at least two accelerometers were mounted 
with one having its sensitive axis oriented parallel to the Z-axis (normal) and the other 
parallel to the Y-axis (lateral). At selected stations, a third accelerometer was oriented 
in the X-axis (longitudinal) direction. On the wings, all accelerometers were alined with 
the Z-axis. Whenever supplementary data were required, additional measimements were 
made with a vacuum -mounted movable accelerometer or by installing additional trans- 
ducers in particular areas and on specific substructures. Strain gages were installed 
on selected skin panels of the fuselage midsection (sidewalls and bottom) and on the 
cargo -bay and cargo -door longerons (see fig. 6) in order to monitor panel stress condi- 
tions during longitudinal preload tests. The gages were motmted on opposite sides of 
the skin and near the center of the panel of interest. Each gage of the back-to-back 
pair was electronically independent of the other. Strain gages were mounted only on the 
outboard flanges of the longerons. 

Data acquisition system .- A schematic of the data acquisition system is shown in 
figure 11. Data signal cables from the transducers were connected to data group switches 
in prearranged combinations so that a group of 14 selected channels of information could 
be monitored or recorded simultaneously. Four data channels were common to all groups: 
the force-gage output, input acceleration (exciter), oscillator (reference) output, and time 
code. The data were also recorded on analog tape for subsequent data reduction. The 
group of instruments designated "On-site analysis" permitted preliminary, quick-look 
evaluation of either real-time or tape-recorded data signals. 

Test procedure .- The experimental data were obtained in basically the same manner 
for all test conditions. The input-force exciter supplied a constant-amplitude sinusoidal 
force to the model at slowly increasing frequency. Selected transducer outputs were 
plotted as diagrams of response amplitude and frequency (unfiltered). These signals 
were also simultaneously displayed as Lissajous patterns of force relative to accelera- 
tion on an oscilloscope, and peak response frequencies and phase shifts were noted. 
Examination of several diagrams of response amplitude relative to frequency served to 
identify individual peak response frequencies, which were then examined by manually 
tuning the frequency of excitation to obtain peak amplitude, dwelling at that frequency, 
and recordii^ the output of all transducers. Apparent modes were also examined more 
closely by using the Kennedy-Pancu method described in reference 23 in order to reveal 
the existence of proximate resonant frequencies and modes. Damping values were cal- 
culated by using the phase-change method outlined by Mead in reference 24. 

Description and results of subsequent static -load deflection measurements of the 
1/8-scale orbiter model are presented in appendix A. These extensive load tests were 
conducted to evaluate the flexibility of various substructures and areas of the model 
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and to provide flexibility influence coefficients for guidii^ revisions to analytical repre- 
sentations which are also discussed subsequently. 

Test Results and Discussion 

The measured natural frequencies in the vertical (±Z) direction of seven responses 
and lateral (±Y) direction of three responses obtained during vibration /tests in the low 
frequency range are listed in table II along with frequencies obtained with the initial 
finite -element analysis (ref. 21). A description of each of the identified modes and the 
dampir^ coefficient obtained for the first mode are included. Initially, wing first bend- 
ing (symmetric mode 2) appeared to be at a different frequency for the left and right 
wir^ components. Subsequent investigation indicated that local responses of wing skin 
panels were masking the response of primary structure and causing each wing component 
to respond differently, due to asymmetries in skin panel conditions. Lightweight angle 
stiffeners were attached diagonally across all rectangular panels (44) between ribs and 
spars to suppress only local skin responses, thereby permitting the identification of the 
normal-mode frequencies which were then found to be nearly coincident. The added 
stiffener weight was 2.4 percent of the original wing weight. 

The resonant frequencies obtained from the initial NASTRAN finite -element model 
were in considerable disagreement with experimental frequencies for all the modes (up 
to 46 percent). Generally, the mode-shape characteristics were similar, as shown by 
comparison of the first four symmetric modes in figures 12, 13, 14, and 15. These 
results indicated that perhaps the experimental model was more flexible than the initial 
analytical representation, that mass distributions were improperly represented analy- 
tically, or that a combination of these existed. There were many possible sources of 
additional flexibility in the physical model such as weak structviral joints, soft connec- 
tions between major substructures or masses (e.g., payload simulation) and supporting 
structure, and model design characteristics such as large -aspect-ratio panels and lack 
of panel stiffening. The investigation of these problems is reported in the section com- 
paring analytical and experimental results. 

Results of vibration tests with tension preloads applied to the fuselage are shown 
in table III. With increase of load to 2670 N (600 lb), there was a small increase in 
first-mode frequency for the experimental model. The results suggested that the fuse- 
lage skin panels, which were observed to have out-of-plane deflections, were not fully 
loaded and hence had reduced stiffness in the model. This possibility was verified by 
strain-gage static data obtained durii^ these preload tests and shown in figures 16 and 17. 
The variation of strain with longitudinal tension load for three typical skin panels is shown 
in figxire 16. As applied load increases, the strain increases on both the inside and out- 
side (opposite) gages of each skin panel. However, the level of strain is different on 
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each side of each panel, which indicates that the skin of the panel is bending and there- 
fore is not fully effective. The differences in strain are quite large except for the 
left-hand sidewall panel. Similar results were obtained from all of the seven panels 
instrumented. 

Figure 17 presents the data obtained from the strain gages on the fuselage cargo- 
bay longerons and the corresponding cargo-door longerons on both sides of the model. 

The differences in strain between door and fuselage longerons throughout the loading 
range indicate that the load paths are not those intended by design. The door longerons 
are not as effective as the fuselage longerons. The door attachment fittings to the fuse- 
lage were suspect in this case. 

Other preliminary test parameter variations were conducted to explore and identify 
suspected flexibility anomalies between the experimental and finite -element models. 

These included the following brief vibration tests: (1) with and without the cargo-bay 
door, (2) with modified experimental boundary conditions at the payload simulation beam 
aft connection to the fuselage, (3) with an increased number of wing -to -fuselage attach- 
ment fasteners (doubled), and (4) with and without the fin ballast. The first three varia- 
tions resulted in little or no change of the first-mode frequency. However, removal of 
the fin ballast increased the first-mode frequency from about 43 to 48 Hz. The amplitude 
of longitudinal motion at various elevations along the aft fin-fuselage surface is plotted 
in figure 18 for the ballast-on and ballast-off conditions. In both cases there is an abrupt 
change in magnitude of motion at the point of fin attachment to the fuselage. Fin motion 
relative to the fuselage is large, which indicated considerable flexibility in this model 
joint. These brief tests indicated that the disagreement in model frequencies between 
experiment and analysis was due to greater flexibility of the physical model. Probable 
sources of reduced stiffness in the physical model were initial out-of-plane imperfections 
in skin panels and several suspected structural joints (tail-fuselage and forward cabin- 
cargo bay junctures were prime suspects). Because of the anticipated cost in time of 
modifying the experimental model and the lack of data needed to positively ascertain the 
major problem areas, as well as the extent of their contribution to the overall correla- 
tion problem, it was decided that the analytical model should be revised to more accu- 
rately represent the existing test structure. The following sections and appendixes 
present discussions and results of both experimental and analytical work conducted to 
achieve better agreement. 


ANALYTICAL MODELS 
John L. Sewall 

The analytical representation of the 1/8-scale orbiter vibration model is a finite- 
element model sufficiently detailed to reproduce overall dynamic behavior. No deliberate 
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effort has been made to simulate local motions, such as panel vibrations. The following 
sections contain brief descriptions of the analytical model, the solution procedure, and 
some experience in its application. The analytical results are compared with experi- 
mental vibration data for the physical model in the concluding section. Further infor- 
mation and related peripheral studies are detailed in appendixes B to E. 

The finite -element analysis used for this model is the NASTRAN Level 15.5 system. 
Development of the analytical model is based on reference 25 and is fully detailed in ref- 
erences 21 and 26. 

NASTRAN Model 

The analytical orbiter model is made up of five substructures as shown in figure 19. 
Fovir substructures - fuselage, wing, door, and fin - appear in figure 19(a) and are each 
represented by combinations of panel, beam, rod, and concentrated-mass elements. The 
fifth substructiire - payload - is simulated by an eight-element beam and concentrated 
masses located along the orbiter center line as shown in figure 19(b). Most of the exter- 
nal surfaces are modeled by membrane elements, which allow in-plane deformations in 
tension, compression, and in-plane shear. Membrane elements are also used for the for- 
ward cabin ballast, forward keel, and the top cover of the wing carry-through structure 
(see fig. 19(a)). Some panels in the cargo-door surface are modeled by plate elements, 
which allow both in-plane and out-of-plane (or bending) deformations. Longitudinal and 
bending stresses are carried by beam elements along the fuselage-door interface. Longi- 
tudinal stresses along the fuselage center line and lower outer edge are carried by rod 
elements. Rod elements are also used, along with shear panels, to model the main keel, 
ribs, and spars in the wing and fin, the aft-engine support structure, the accordionlike 
door frames, and the fuselage cross-sectional frames. The door-frame simulation is 
shown in figure 19(a) and is intended to prevent transmission of longitudinal stresses 
through the door surface. Typical fuselage frames are shown in figure 20 for the cabin 
area, payload bay, and engine -support area. 

Solution Strategy 

The analytical natural frequencies and mode shapes were calculated by solution of 
the general eigenvalue equation 



( 1 ) 


where is the stiffness matrix, [j\^ is the mass matrix, X is the eigenvalue, 
and {u^ is the displacement vector, from which the eigenvectors are obtained. The 
mode shape (or eigenvector) is a set of relative, or normalized, displacements. 
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In this investigation, the mode shape was normalized to the maximum relative dynamic 
displacement in the structure (i.e., <l>i = Uj^Uj^^x every point i of a given mode]. 
The sizes of and are governed by the number of unconstrained degrees of 

freedom at grid points in the structtire. 

The large number of degrees of freedom associated with the grid selected to repre- 
sent the entire orbiter -model structure precluded a solution of equation (1) in a sii^le 
NASTRAN computer execution on the NASA Langley Research Center CDC computer. 
Instead, the complete analytical model was divided into five major substructures as 
shown in figure 19. An adaptation of the current operational substructuring procedure 
was used, as described in reference 26. The adapted substructuring procedure detailed 
in reference 26 is contained in a single computer -program alteration in each phase of 
the calculation. 

As indicated in the simplified block diagram of figure 21, there are four sequential 
operations involved. In phase I of the calculation, stiffness and mass matrices are gen- 
erated for each substructure and stored on tapes. This information is then recopied 
onto a single tape in the so-called tape-copy operation and is used in phase II to generate 
the natural frequencies and mode shapes of the assembled orbiter. Phase II is designated 
"pseudo-structure" because its grid network includes only enough of the complete orbiter 
network to properly define structural dynamic behavior of the interfaces between sub- 
structures and overall orbiter mode shapes. More detailed mode-shape distributions 
could be obtained in phase in of the calculation, but this is an optional operation and was 
not executed in this investigation. Substructuring is discussed further in appendix B and 
in reference 21, vol. II, and reference 26, and the relation of substructuring to analytical 
modeling philosophy is also discussed in references 21 and 26. 

Operational Experience 

As is evident in figure 21, a sir^le calculation of orbiter frequencies and mode 
shapes by the substructuring procedure requires a total of seven separate computer 
executions (or runs), not including phase III. A measure of the size of each run is 
given by the typical central processing unit (CPU) time indicated in the figure for carry- 
ing out the eigenvalue analysis on the NASA Langley Research Center CDC-6600 series 
computers. The phase I fuselage and phase II orbiter runs used the most CPU time and 
required CDC 6600 machines with nearly all the available storage capacity (i.e., field 
length). A more complete listing of computer requirements is given in table IV. The 
phase II runs required a field length of 300 000 octal words. The other runs required 
less time and field length, and the smallest ones (tape copy and Phase I payload) could 
be processed on the CDC 6400 with small amounts of CPU time. The number of runs 
gives an indication of the scope of the investigation into factors affecting the correlation 
between analysis and experiment. More variations in parameters were made for the 
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analytical symmetric model than for the antisymmetric model. The cumulative CPU time 
of seven runs, resulting in an orbiter phase II run, was 7915 sec, as shown in figxire 21. 

Another measure of size in the calculated procedxire is evident in the numbers of 
elements, grid points, and degrees of freedom listed in table V, which is taken from 
table 2 of reference 21, vol. I. The most numerous elements are rods, membranes, 
and shear panels - in that order. The maximum matrix size in equation (1) for any 
substructure or for the entire orbiter may be estimated by simply multiplyir^ the num- 
ber of grid points by the maximum number of degrees of freedom possible at each grid; 
this maximum number is six, three displacements and three rotations. For example, 
the phase I fuselage can have at most 3222 degrees of freedom, and the phase II orbiter 
can have a maximum of 5604 degrees of freedom. However, the imposition of constraints 
results in a reduced number of degrees of freedom. In NASTRAN, the constraints are of 
two kinds: (1) single point constraints (SPC) restrictii^ certain degrees of freedom at 
grid points so as to satisfy boundary conditions and symmetric or antisymmetric condi- 
tions, and (2) multipoint constraints (MPC) relating motions between different parts of 
the structure, thereby altering the number of independent degrees of freedom. The 
numbers of degrees of freedom remainir^ in the present investigation after imposing 
SPC and MPC are so indicated in table V. Because of model symmetry, certain SPC 
restrict motions to vertical, longitudinal, and pitch in the vertical plane of symmetry 
throxigh the orbiter center line, and other SPC allow only lateral, roll, and yaw motions 
along this center line; hence, the two models, symmetric and antisymmetric. For the 
phase I fuselage and Phase II orbiter, the reduction in maximum matrix sizes is well 
over half. This is also true for the phase I wing and door, but for the phase I fin and 
payload, the reductions were about three to one. i 

Table V also indicates another substantial reduction in the number of degrees of 
freedom through application of the Guyan reduction (ref. 27). In the present investiga- 
tion, selected degrees of freedom were reduced out by this procedure from all substruc- 
tures except the phase I payload. For phase n, the reduction in degrees of freedom was 
not so great as for the individual substructures but did result in dynamic matrices of 
order 339. This problem required from 3900 to over 4000 sec of CPU time to obtain 
nine natxiral frequencies and mode shapes by the inverse power method on the Langley 
Research Center CDC 6600. 

Without the substructuring procedure, the NASTRAN solution of equation (1) for 
the orbiter was intractable. With substructuring, the limits of present computer cap- 
ability were closely approached for practicable operation in the present investigation. 
Priorities at the Langley computer complex were essential for reasonable turn-around 
times (e.g., overnight). Further aspects of the operational experience with the present 
l/8-scale dynamic model are discussed in appendix B and in reference 21, vol. I. 
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static Analysis 


In conjunction with the static test program described in appendix A for the physical 
model, a NASTRAN analytical model was formulated to compute displacements due to 
static loading from the relation 



where is a vector of static discrete loads Pj located in various parts of the model. 

Solutions for the displacements {“) were obtained with Pj = 4.45 N (1 lb). These cal- 
culations were performed in a single computer execution having as input a simple com- 
bination of all five substructures and involvir^ no substructurirg operations. Over 
2400 degrees of freedom were involved, as indicated by the totals for the columns enti- 
tled "After SPG and MPC" in table V. Each of these calculations took about 2800 sec 
of CPU time and were run at 300 000 octal words of core. In an effort to gain as much 
as possible in computational efficiency, the NASTRAN grid point resequencing feature 
was used to give well-banded matrices (i.e., matrices each with nonzero elements clus- 
tered about the main diagonal). The SEQGP data cards required for this operation were 
generated automatically by means of the BANDIT computer program described in refer- 
ence 28. Results checked satisfactorily with static deflections calculated at GAG using 
the orbiter phase II pseudostructure (ref. 21, vol. II). Other aspects of the operational 
experience with the 1/8-scale static model are discussed in appendix B and reference 21, 
vol. II. 


ANALYTIGAL RESULTS AND GOMPARISON WITH EXPERIMENT 

John L. Sewall 

Unacceptable discrepancies revealed in initial comparisons between analytical and 
measured frequencies (table II) prompted an investigation resulting in the identification 
of significant structural joint and panel behavior not accounted for in the initial mathe- 
matical modeling. The investigation involved studies of the effects of various modifica- 
tions to the NASTRAN model: first, as reported herein, second, in supporting studies 
in appendixes A, G, D, and E, and third, in reference 21. Primary emphasis is given 
to symmetric model conditions, although some results are included for the antisymmetric 
model. These studies led to various revisions in the analytical model, and the frequen- 
cies of the best of these models are compared in table VI with measured frequencies. 
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Initial and Side Investigations 


The initial and side investigations are listed as follows (where the tabulation deals 
with items which were found to be of minor importance to the analytical-experimental 
correlation): 

Fuselage Satisfactory agreement between 

NASTRAN and section bending 
stiffness (ref. 21, vol II) 

Generalized mass estimates Fuselage and ballasts are largest 

contributors (ref. 21, vol. II) 

Cabin ballast modeling Negligible effect of replacir^ mem- 

brane by plate elements 

Wing carry -through modifications Negligible effects of reducing 

panel thickness to simulate cut- 
outs and shear connections 

Door longeron effectiveness Negligible dynamic effects of 

reduced door longeron cross- 
sectional area to account for more 
bendii^ flexibility (appendix C) 

Fin ballast support Uncertain effects of local distortions 

Panel modeling accuracy Satisfactory accuracy of NASTRAN 

membrane panels (appendix D) 

High panel aspect ratio Could cause as much as 5 -percent 

increase in frequency (appendix D) 
and ref. 21) 

None of the items in the foregoing listing revealed the major cause of the 
experimental -analytical frequency discrepancy, although they did lead to a better under- 
standing of both the physical and analytical models. Further elaboration of these items 
follows. 

Fuselage .- Satisfactory agreement was obtained for bending stiffnesses extracted 
from the NASTRAN generated stiffness matrix compared with bending stiffnesses deduced 
from a NASTRAN generated deflection curve and calculated from section geometry of the 
physical model. 

Generalized mass estimates .- Generalized mass estimates based on the first sym- 
metric mode helped to identify areas of the model requiring further investigation. 
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Dominant portions of the generalized mass existed in the fin ballast, the orbitii^ propul- 
sion system (OPS) ballast, cabin ballast, the remainder of the fuselage, and the wing. 

Cabin ballast remodelii^ . - Negligible reductions in analytical frequencies were 
obtained by allowing more degrees of freedom in the cabin ballast and by attempting to 
simulate effective stiffnesses in the wing carry -through structure and door longeron. 

The ballast was in the form of a thick (0.95 cm (3/8 in.)) plate attached to the cabin 
structure, and its presence could conceivably change stiffnesses in this area. Replacii^ 
membrane elements with plate elements in the cabin ballast simply allowed additional 
pitch and roll degrees of freedom not permitted with the membranes. 

Wing carry-through modifications .- Views of the wing carry -through structure 
including pertinent interface details are shown in figure 22. Figure 22(a) is an overall 
view of the four main wing -fuselage connections. Figure 22(b) is a closeup of the con- 
nection at fuselage station 144.75 and shows the cutout in the top-cover panel for the 
fuselage frame and also the vertical channel member to which the wing -fuselage shear 
pin connection is attached. Membrane panel thicknesses in the top cover of the 
NASTRAN model (fig. 19(a)) adjacent to the wing -fuselage interface were reduced to 
account for a possible loss in wing bending stiffness due to the cutouts. In the initial 
NASTRAN model, the thickness of each vertical shear panel in the wing root shear con- 
nection was arbitrarily chosen to be the same as the channel thickness (fig. 22(b)). In 
the present investigation, this thickness was reduced according to the cross-sectional 
area of the shear pin to approximate a more realistic shear flow between wing and 
fuselage. 

Door longeron effectiveness .- A typical door-to-door fuselage longeron attachment 
is the door clip shown in figure 6, and there are four of these nearly equally spaced along 
the cargo bay. Longitudinal loads transmitted by shear pins through the clips from the 
fuselage longeron to the door longeron will result in a combined longitudinal and rota- 
tional (pitch) deformation of the door longeron. However, the bending flexibility in the 
door longeron due to this rotation had not been considered in the initial NASTRAN model. 
When it was considered, as shown in appendix C, the cross-sectional area of the door 
longeron was approximately halved to account for a reduced door longeron effectiveness. 

Fin ballast support . - During vibration tests of the physical model, some local 
deflections of the fin-ballast support structure could be deduced from the vibration mode 
shape shown in figure 18. Attempts to reproduce these local motions by revisions to 
the analytical model involved refinements to the fin-ballast support structure as illus- 
trated in sketches (a) and (b). These refinements consisted of a NASTRAN model of 
the channel-shaped support structure and alternative ways of connecting it to the fin 
ballast. In the initial model, the motions of the fin ballast were simply linked by MFC 
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relations to the motions of the grid network of the fin-tip rib. Sketch (a) indicates one 
model revision in which the MPC links are confined to two points along the midplane of 
the fin. Two variations of this model resulted in drastic reductions in the phase I fin 
frequencies and imrealistic mode-shape distortions along the fin-tip midplane. The other 
revised model in sketch (b) allowed for additional MPC links, as indicated, to three points 
along the channel flange but resulted in negligible phase I frequency reductions and mode- 
shape changes compared with those of the initial model. Consequently, modifications to 
the initial ballast support structure model were considered unwarranted. 

Panel modeling accuracy .- The capability of NASTRAN to predict accurately the 
dynamic response of a fuselagelike structure was evaluated in a side investigation 
reported in appendix D. Frequencies, mode shapes, and dynamic shear stresses of 
NASTRAN analyses of a simplified model of the cargo bay agree well with those obtained 
in two closed-form solutions. 

High panel aspect ratio . - As is evident in figure 19(a), many of the analytical mem- 
brane panels in the fuselage sidewalls and bottom are long and narrow and, therefore, 
possibly have too much in-plane stiffness. This characteristic is shown in reference 29 
for in-plane load-displacement studies of NASTRAN quadrilateral plate elements. As 
shown in reference 21, vol. II, and in appendix D, this high-aspect-ratio effect is esti- 
mated to account for an approximately 10-percent increase in stiffness and thus an 
approximate 5 -percent increase infrequency for NASTRAN quadrilateral membrane 
elements (QDMEM2). 


Effects of Joints 

The sensitivity of the model’s vibration behavior to joint flexibility is indicated as 


follows: 

Fin-fuselage interface Accounted for about 13^ percent of frequency 

discrepancy (appendix C) 

Cabin-cargo bay Negligible effect on frequency 

Payload supports Small effect on lower mode frequencies 

(ref. 21, vol. II) 
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As shown in appendix C, joint flexibility estimates were determined both by structural 
analysis of the actual connections in the physical model and by correlation of measured 
and calculated static deflections. 

Fin-fuselage interface .- As shown in figure 23, the fin is attached to the fuselage 
at the front, center, and rear spars of the fin and by angle splices along the fin-fuselage 
junction. As previously described, the fin fore and aft spars are connected by metal clips 
(visible in fig. 23) to fuselage frames at the ends of the engine support structure, and the 
center fin spar is connected by a clip and a gusset (fig. 23(a)) to a fuselage frame in the 
middle of the engine support structure. In between these three fuselage frames, the fin- 
fuselage angle connection is fastened only to the skin without any underlying backup or 
support structure. 

NASTRAN models of the fin-fuselage interface are shown in figure 24. In all cases, 
the actual connections were modeled by grid points along the fin center line and the out- 
side edge, and the motions of both substructures were matched at these points. There 
were no springs in the initial NASTRAN model. Phase II orbiter initial analyses at GAC 
showed a 5 -Hz frequency change due to alterii^ boundary conditions of the forward fin- 
fuselage attachment and a negligible change due to altering the aft fin attachment (see 
ref. 26). In addition, deflections due to static loadii^ measvired along the fin trailing 
edge and aft end of the fuselage revealed considerable flexibility in the overall fin-root 
attachment not accounted for in the original NASTRAN model. This flexibility was 
modeled in different ways in terms of linear spring constants, as illustrated in figure 24. 
The 13^-percent frequency discrepancy noted in the second of the foregoing tabulations 
applies to the first symmetric mode which was calculated with flexibility allowed in the 
forward fin attachments only (model F2, fig. 24) and based on the difference between 
measured and analytical static deflections, as shown in appendix C. 

Cabin-cargo bay and payload supports .- Joint flexibilities were also investigated 
and simulated at the cabin-to-cargo bay juncture (appendix C), shown schematically in 
figure 25, and at the payload mountir^ points in the fuselage (see ref. 21, vol. II). Neither 
of these joints had as large an effect on fundamental frequency as the fin-root flexibility. 

Panel Effectiveness 

Panel load-carrying effectiveness is listed as follows (where the tabulation iden- 
tifies the two types of modeling chaises that accounted for the major part of the 
experimental-analytical frequency discrepancy, namely, partial effectiveness of exte- 
rior surface panels in transmitting inplane direct and shear stresses): 
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Rod-shear panel naodifications Accounted for 53 percent to 80 percent of 

frequency discrepancy (appendix E and 
ref. 21, vol. II) 

Reduced membrane -panel thicknesses .... Accounted for about 80 percent of frequency 

discrepancy (table VI) 

In the first type of panel model, membrane elements were replaced by shear panels 
and bordered by rods. The second type of effective panel model was simply the membrane 
element with its thickness reduced according to criteria derived from independent studies 
described both in appendix E and reference 21. A thickness reduction in the membrane 
panel reduces direct and shear -stress effectiveness simultaneously, whereas these prop- 
erties can be separately varied in the rod-shear panel model. 

Rod-shear panel modifications .- As shown in appendix E, rod-shear panel modifi- 
cations were made for the cargo-bay part of the fuselage and for top and bottom covers 
of the wii^ and wing carry-through structure. Sizing of rod cross-sectional areas to 
be effective in direct stress was based on correlation of orbiter measured and analytical 
static deflections, and the choice of shear -panel thickness was made independently, and 
arbitrarily, to represent partial effectiveness of the panels in shear. Six different 
models were investigated and are designated RSI to RS6. The best rod-shear panel 
model is RS6, which is model II of reference 21. Results from this model are com- 
pared with those of the other rod-shear panel models in appendix E. 

Reduced membrane panel thicknesses .- Reductions in membrane panel thicknesses 
were based on results of independent studies that included theoretical criteria based on 
reference 30 and tests and finite -element analyses of a simplified cantilever panel model 
approximately the same size as the cargo-bay sidewall (see also appendix E). In con- 
trast to the rod-shear panel models, thickness reductions were made for all exterior 
surface membrane panels of the fuselage, wii^s, and the top cover panel of the wing 
carry -through structure. Further discussion of these modifications is also given in the 
following sections and in appendix E. 

Symmetric Natural Frequencies and Mode Shapes 

Symmetric natural frequencies showing the combined effects of joint flexibility 
and reduced panel effectiveness are listed in table VI(a) for the best analytical models. 
Measured frequencies are also included for the first four modes. Mode identifications 
are given in table II. 

For the rod-shear panel model RS6, magnitudes of the reduced effectiveness in 
direct stress were obtained by application of the large -deflection equations for panel 
imperfections in reference 30 and amounted to 50 percent for the fuselage payload-bay 
region and 85 percent for the wing. The payload-bay panels were considered to be 
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two-thirds effective in transmitting shear and the wing panels were 85 -percent effective. 
These values were chosen somewhat arbitrarily, and the choice of different values for 
fuselage and wing was influenced partly by different panel aspect ratios in wing and fuse- 
lage. Different joint flexibilities were used in reference 21 than in the models considered 
in the present paper, and additional flexibilities were introduced at the payload supports. 

In the reduced membrane model RMl, a panel effectiveness of 60 percent was 
uniformly applied in both direct stress and shear to all exterior membranes of the fuse- 
lage and wings, including the wing carry -through cover panels. This value was based 
partly on the criteria used for model RS6 and partly on results of some simplified panel 
studies described in appendix E. Frequencies of the reduced membrane model are in 
close agreement with those of the rod-shear panel model RS6, and both sets of frequen- 
cies generally agree better with measvired frequencies than do those of the other rod- 
shear panel models considered in this study. 

Analytical symmetric mode shapes for model RMl are compared with measured 
mode shapes in figures 26 to 29. In general, the agreement is about as good as that 
for the initial NASTRAN model in figures 12 to 15. Mode shapes of the wing and fuse- 
lage in the first mode are in somewhat better agreement for RMl than for the initial 
model, but just the reverse is true for the payload and fin. Better agreement with 
experiment is also shown for RMl in the outer wing for the third mode and in the fin 
for the fourth mode. 

Antisymmetric Natural Frequencies and Mode Shapes 

Analytical frequencies of the best antisymmetric models are included in table VI(b). 
Corresponding mode shapes are shown, together with experimental mode shapes, in fig- 
ures 30 to 33. In addition to the 60-percent panel effectiveness in model RMl, joint 
flexibility is included for both the cabin-to -payload bay junction and for the fin-fusel^e 
interface using a roll-spring model (F3 in fig. 24) similar to one derived in reference 21. 
The frequencies of the model are lower than those of the initial antisymmetric model by 
about the same amounts as those of the symmetric model; and the first antisymmetric 
frequency of model RMl is in as good agreement with the measured first antisymmetric- 
mode frequency as are the symmetric -mode frequencies. 

The fourth antisymmetric frequency (77.8 Hz in table VI(b)) is associated with a 
predominant roll (or torsional) response of the orbiter. Much of the torsional stiffness 
of this mode can be attributed to the closed doors of the orbiter. The presence of a low- 
frequency torsional mode in the phase I fuselage substructure listed in table Vn(a) sv^- 
gests the possibility of a low orbiter (phase II) torsion mode with the doors open, and this 
possibility is confirmed both by a phase II calculation with the doors off and by test, as 
shown in table Vll(b) and figure 33. Modes 2 to 6 in table VH(b) correspond to modes 1 
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to 5 in table VI(b), and their frequencies are increased from 2.3 percent to 14.8 percent 
due to loss of door mass. A low-frequency torsion mode was also obtained for the highly 
simplified cargo -bay model of appendix D. 

Experimental antisymmetric mode shapes with doors on are compared in figures 30 
to 32 with correspondii^ analytical mode shapes for the first, second, and fifth modes. 

t 

The agreement is generally good except for the large discrepancies evident in figiires 30 
and 32 for the lateral payload component. This disagreement may be due to inadequate 
analytical representations of payload stiffnesses and lateral flexibilities in the physical 
constraints at fuselage -payload interfaces. The boxlike structure of the payload beam 
is formed from two channels intermittently welded at the edges of their flanges in a 
series of discontinuous joints, each 5.34 cm (2.1 in.) in length, along the beam, with 
additional stiffness provided by plates attached to the vertical sides (or channel webs) 
of the box. Lateral bending and torsional stiffnesses of this beam for the analytical 
model are based on the assumption of continuous welds and may be too high. However, 
a realistic approximation of these stiffnesses is not readily available. Moreover, arbi- 
trary reductions in stiffness have shown small reductions in phase II frequencies and 
negligible changes in the analytical payload lateral mode shape. 

Another relatively minor mode-shape disagreement in the antisymmetric vibration 
modes exists along the door longeron. Here, in the range of about 93 to 105 Hz, the door 
longeron experienced numerous large lateral displacements that not only did not match 
the smoother overall deflections of the adjacent fuselage longeron, but also could not be 
reproduced in the present analysis. A closer analytical representation of these motions 
would require a finer grid network for the door, together with the allowance of a suffi- 
cient number of lateral degrees of freedom along the door longeron. 

CONCLUDING REMARKS 

Vibration tests and analysis of a 1/8-scale model of a preliminary design for a 
space shuttle orbiter have been reported. Test results for four symmetric vibration 
modes and four antisymmetric modes have been compared with finite -element predic- 
tions based on NASTRAN analyses. Initial analysis results using geometrically derived 
stiffnesses indicated poor agreement with test data. 

An extensive analytical and experimental investigation of the structural response 
of the orbiter, consistii^ of parallel statics and vibration studies, revealed that the 
causes for the discrepancy were (1) out-of-plane initial imperfections of the panels in 
the fuselage (about 85 -percent contribution) and (2) a soft connection between the model 
fin and fuselage (about 15 -percent contribution). These characteristics of the physical 
model could not be defined prior to its construction. When these characteristics were 
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included in the mathematical model, satisfactory correlation was obtained; symmetric 
and antisymmetric frequencies could be predicted within 11 percent for the modes 
compared. 

The conclusions from parallel studies contributir^ to this good correlation should 
be noted. The static test-analysis program was essential in determinii^ realistic esti- 
mates of flexibility in various parts of the present structure where initial analytical pre- 
dictions were overly stiff, such as in the fin-fuselage joint and the cargo-bay doors. 
Comparisons of analytical results for coarse and'fine grid models of a simplified orbiter 
fuselage with closed-form solutions confirmed the NASTRAN accuracy in predicting nat- 
ural frequencies and mode shapes of monocoque structures. This study also showed that 
the use of relatively high aspect ratio of the NASTRAN membrane elements does not 
result in a frequency increase of more than a few percent for the present structure. Ana- 
lytical and experimental panel-effectiveness studies indicate that panels such as those 
forming the fuselage and wing skins are about 60-percent effective in carrying direct 
stress and in-plane shear over an aspect-ratio range of 0.5 to 2 and having initial out-of- 
plane imperfections of one fourth the skin thickness or greater. 

The NASTRAN finite -element program is capable of satisfactorily predicting 
dynamic characteristics of structures similar to the space shuttle orbiter. However, 
even with the use of the sub structuring capability of the program (in Level 15) and pitch - 
plane symmetry, a large structural dynamic mathematical model resulted which utilized 
most of the LRC computer storage capability and required long run times. Two and one- 
half hours were required on the CDC 6600 to obtain the lowest symmetric and antisym- 
metric modes (four each) after reduction of the assembled model (2500 degrees of free- 
dom) by substructuring and Guyan elimination to 332 degrees of freedom. Thus, 
parametric studies may become time consuming and expensive. 

Langley Research Center 

National Aeronautics and Space Administration 
Hampton, Va. 23665 
July 21, 1975 
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APPENDIX A 


STATIC -LOAD DEFLECTION MEASUREMENTS 
Ulysse J. Blanchard and John E. Flynn* 

Introduction 

Initial experimental and analytical dynamic response data obtained for the l/8-scale 
shuttle orbiter model showed a significant discrepancy in natural frequencies. Values 
obtained by analysis were higher than by experiment. The data also indicated that flexibil- 
ity of the experimental model was greater than the originally formulated NASTRAN finite - 
element representation. It appeared that a combination of partially effective skin panels, 
due to deviations in flatness, and joint flexibilities could be contributors to the disagree- 
ment. Extensive static -load tests were conducted in order to evaluate the flexibility of 
various substructures and areas of the model and to provide flexibility influence coeffi- 
cients for guiding revisions to analytical representations. 

Apparatus and Procedure 

The static-load tests were conducted with the model supported at interstage connec- 
tion fittings. The photograph of figure 34 shows the model mounted on four pedestals 
which were attached to steel beams embedded in the concrete floor of the Langley struc- 
tural dynamics research laboratory. Photographs of the orbiter interstage fittings and 
adapter brackets attached to the pedestals are shown in figures 35(a), 35(b), and 35(c). 
Support point constraints of the orbiter during static -load testing are illustrated in fig- 
ure 36 and these essentially duplicate conditions for the orbiter mated to the external 
tank. 

The points of application and the direction in which symmetric and antisymmetric 
loads were applied to the model are illustrated in figures 37(a) and 37(b). In all, 14 sym- 
metric and 4 antisymmetric load conditions were investigated. Some of these conditions 
also were investigated with and without the cargo-bay doors. Incremental loads were 
applied with calibrated weights placed on various weight pans connected through cables 
to primary structural members. Some of these loading apparatus are shown in fig- 
ures 34 and 38. 

Dial indicators on support frames (figs. 38(a) and 38(b)) were used to measure 
deflections of the model resulting from the applied loads. The basic dial -indicator 
coverage of the model durir^ the tests is illustrated in figure 39. Measurements were 

*Grumman Aerospace Corporation. 
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made at points on the model having maximum stiffness such as the fuselage frames, 
longerons, juncture of wing spars and ribs, and ballast masses. 

During the tests, several loading and readout procedures were employed to mini- 
mize error and to maintain acceptable data accuracy. As an example, consistency in 
the data was improved by cycling an intermediate load on and off the model several times 
prior to each test run. Early in the investigation a statistical check of measurement 
accuracy was made. Results are shown in table VIII. Nineteen on-off loads simultane- 
ously applied at two weight pans at the nose of the model resulted in a 3 -sigma deviation 
of about 12^ percent from mean for deflection readings. 

Results and Discussion 

Typical load deflection data obtained during loading of the aft fuselage structure 
are shown in figure 40. Measurements obtained at the loading point, the fin tip, and the 
lower fuselage corner are shown for longitudinal loads applied to ballast masses repre- 
senting the orbiting propulsion system (OPS). Data plotted above the zero -load line are 
for loads applied in the aft (+x) direction and the data below the line are for loads applied 
in the forward (-x) direction. As shown in the figures, straight line fairings of the data 
points were made whenever possible. Generally, the loading and unloading legs for each 
test run had to be treated separately due to various degrees of hysteresis. Also, the 
individual loading directions were faired separately since in most cases there was some 
nonlinearity going through zero load. The best data were consistently obtained in areas 
of maximum deflection such as the loading point. As distance of the measurement point 
from the load point increased or as total deflections decreased, the measurements became 
more erratic. 

Slopes of the faired lines were determined and then averaged for each load appli- 
cation point. The values obtained are presented in table IX as measured flexibility coeffi- 
cients at the load points. Also listed in the table for comparison are deflection coeffi- 
cients calculated using the initial NASTRAN analysis (model I of ref. 21). 

Concludii^ Remarks 

In all test cases shown, the experimental model was more flexible than the original 
finite -element representation. The data for symmetric loading at midfuselage (bending) 
show that the experimental-model cargo door was less effective than predicted by analy- 
sis. This complements strain-gage data obtained during previous vibration tests (fig. 17). 
Also, a much larger flexibility in the experimental fin to fuselage joint is indicated by 
results for symmetric loading of the fin ballast. Another significant result is the greater - 
than -predicted torsional effectiveness of the experimental cargo door (antisymmetric nose 
loading case). 
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NASTRAN OPERATIONAL EXPERIENCE 
Barbara J. Dvirlir^ 

This appendix presents various aspects of the operational experience with the 
CDC version of NASTRAN Level 15.5 during the dynamic and static analyses of the 
1/8- scale shuttle orbiter model. The discussion records some of the user options 
selected and problems encountered and circumvented. In addition, several comments 
and suggestions are included to benefit the new or occasional user or programer of 
NASTRAN. 


Substructuring 

The analytical natural frequencies and mode shapes for the orbiter were obtained 
by coupling the five orbiter substructures in a substructure analysis procedure in the 
Normal Mode Analysis (Rigid Format 3) in NASTRAN. This procedure utilized an 
adapted DMAP (Direct Matrix Abstraction Program) alter packet in Rigid Format 3 and 
a DMAP tape copy program which were provided by the Grumman Aerospace Corporation 
(see refs. 21 and 26). The single alter packet is easier to use than the NASTRAN sub- 
structure analysis procedure described in section 1.10 of reference 25 which requires 
a separate alter packet for each phase of the substructuring process. Provisions for 
makir^ equilibrium checks and for calculating substructure frequencies and mode shapes 
in phase I, either for the free-free boundary condition or with the substructure fixed at 
an interface and free elsewhere, are also included in this sii^le alter packet. These 
additional features were particularly useful in the orbiter study to assess model changes. 

The DMAP tape copy program was developed for this investigation under the rules 
of NASTRAN's DMAP (described in ref. 25) for creating special and analytical programs 
or routines within the NASTRAN system. This DMAP program is described in refer- 
ence 26 and is executed in NASTRAN after the phase I executions have been completed 
(see fig. 21). It reads the stiffness and mass matrices from each of the tapes generated 
in the phase I runs and places these matrices on a single tape which will be read in the 
phase II execution. 

During this analysis a coding error was discovered in the CDC version of NASTRAN 
Level 15.5 which precluded use of the INP9 user tape. This error was circumvented by 
changing INP9 to INP8 in the DMAP alter packet. 
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GINO Buffer Size 

The GINO (General Input/Output) bxiffer is the storage reserved in open core for 
each GINO file opened. The default GINO buffer size in the CDC version of NASTRAN 
Level 15.5 is 666. When an increased buffer size is desired for CDC NASTRAN runs, it 
is recommended that multiples of 512 be added to the default buffer size. The increased 
buffer size used on some of the orbiter runs was 2202. It was obtained by insertir^ a 
NASTRAN card 

NASTRAN BUFFSIZE = 2202, NLINES = 35 

before the ID card (identification card which is the first card of the NASTRAN Executive 
deck) in the NASTRAN Data deck. The print control parameter modification was included 
so that only 35 lines of data would be printed on each page of 21.6 -cm ^8~in.^ tabulatii^ 
paper. 

Buffer size on substructuring runs .- If an increased buffer size is required on any 
run in the sequence of runs for coupling the substructures, then all runs in the sequence 
must use the increased buffer size. This point is emphasized so that a new user will 
consider the need for an increased buffer size on phase II of the substructuring sequence 
before submitting phase I runs. (See fig. 21 and the section entitled "Solution Strategy" 
in the main body of this paper.) 

Bviffer size on the restarts . - It is mandatory that the GINO buffer size used on an 
initial checkpointed run be maintained on a subsequent run which utilizes as an Old 
Problem Tape (OPTP) any reel generated as a New Problem Tape (NPTP) on any run 
in the sequence of runs. Thus, either the default or modified GINO buffer size can be 
changed only by submitting an initial run for the configuration. 

Working Core 

NASTRAN efficiency should be maximized whenever possible. To achieve this end, 
sufficient working core must be allocated or excess core eliminated. In many applica- 
tions, modules that perform symmetric matrix decompositions, such as RBMG2 and SMPl, 
require the largest amount of core. The core requirements for symmetric matrix decom- 
position are described in section 3.5.14 of reference 31 and also in section 2.2 of refer- 
ence 32. The working core W required for single precision symmetric decomposition 
is given by equation (9) on page 2.2-4 in reference 32 as 


o2 ^2 

WgBR-i|-+BC + ^ 


B g R 5 2 


(Bl) 
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where B is the bandwidth, R is the number of columns of terms inside the band that 
will fit in core, and C is the number of active columns. 

For efficient operation, sufficient core should be allocated so that a spill sitixa- 
tion is avoided. The no-spill situation for symmetric decomposition is indicated from 
NASTRAN user information message 3023 (which prints values of B, C, and R) 
when R = B - 1. The working core W, as given in equation (Bl), does not include the 
storage required for GINO buffers. As previously mentioned, space sufficient for three 
buffers is required when no spill occurs and space sufficient for five buffers is needed 
when spill occurs. 


Diagnostics 

Optional diagnostic output may be obtained by inserting DIAG cards in the NASTRAN 
Executive deck. The DIAG's used most frequently on the orbiter runs were 1, 8, 11, 13, 

14, 15, 16, 19, 21, and 22. A description of the output triggered by these DIAG's is given 
in reference 25, page 2.2-4. Since the reason for selecting some of these DIAG's may 
not be apparent to the new NASTRAN user, brief comments will be given. 

Open core length is printed in the dayfile (program chronology at the end of a com- 
puter printout on the CDC computers) when DIAG 13 is included. This information, used 
in conjunction with the location for the beginning of open core (see ref. 31, pp. 5.5-17 
to 5.5-47), enables the user to determine more accurately the field length actually needed 
for the job. 

Although included primarily so the systems programer can trace GINO OPEN/CLOSE 
operations, DIAG 15 can also serve the user. A significantly large number of these 
OPEN/CLOSE messages in the dayfile indicates that a larger field and/or an increased 
buffer size is needed for the job. When solving very large problems, these messages 
may be an indicator that solution of the given problem either is not feasible or perhaps, 
in the case of a job abort, is impossible due to excessively long run times or insufficient 
core on a given computer. It may be necessary to remodel part or all of the structure 
in order to reduce the number of degrees of freedom allowed and still maintain adequate 
representation of the structtire. 

Run times can be estimated more accurately for subsequent submittals by checkii^ 
execution times for modules. These times can be computed from the BEGIN and END 
times printed in the dayfile for each functional module. Tabulations and plots of these 
run times against problem size for selected modules were useful in the present analysis. 

Grid-Point Resequencing 

A large percentage of the total computing time for structural analysis in NASTRAN 
is associated with triangular decompositions. The decomposition routines treat all 
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matrices as partially banded, that is, nonzero terms clustered near the diagonal are 
treated inside a band of constant width, and nonzero terms outside the band are treated 
separately and referred to as active columns. Since the structural matrices formed in 
NASTRAN are both symmetric and sparse, computational efficiency can be obtained best 
when using the displacement method by selectii^ the numbers assigned to grid points in 
such a manner that stiffness matrices will be created with relatively narrow bands. (See 
ref. 25, section 1.2.2. and ref. 32, section 2.2.) 

External grid-point numbers may be selected in a convenient manner and then a 
NASTRAN option may be used to resequence the external grid points internally. This 
option is activated by the presence of a set of SEQGP cards in the Bulk Data deck. These 
cards provide the correspondence between the external grid-point numbers for the model 
and the internal grid-point numbers that will be used in the calculations. 

Since selection of the internal grid-point numbers for optimizing the bandwidth can 
be difficult and time consuming, programs for automatically generating the SEQGP cards 
have been developed. The BANDIT computer program described in reference 28 was 
used to generate the SEQGP cards used in the static analysis of the orbiter model. This 
orbiter model contained all five substructures, but the substructuring procedure was not 
utilized in the Static Analysis (Rigid Format 1) executions in NASTRAN Level 15.5. 

When the Guyan reduction (ref. 27) featare is being utilized in a Normal Mode 
Analysis run (Rigid Format 3) for a large problem, relatively inexpensive computer runs 
can be made to determine if resequencing reduces the bandwidth. One of these runs may 
involve no re sequencing. Each of the runs should be checkpointed. Each run should be 
submitted with sxafficient time to allow entry in the SMPl module but insufficient time for 
completion of this module. The SMPl module partitions the constrained stiffness matrix, 
solves for the transformation matrix, GO, and performs a matrix reduction to obtain the 
reduced stiffness matrix, KAA (see ref. 25, p. 3.4-7). User information messages are 
output with values for B, C, and R identified in equation (Bl) and the decomposition 
time estimate in seconds. These messages are followed by a system fatal message 
which states that there is insufficient time remaining for RSPSDCM and gives the time 
estimate in seconds (RSPSDCM is the routine used for a real, single precision symmetric 
decomposition in the CDC version of NASTRAN). Following output of this message, the 
job is terminated by NASTRAN. The calculations are then continued on an immodified 
restart run from the run selected. The computer time lost is that spent on runs which 
are not continued. 


Eigenvalue Extraction 

The need for faster eigensolution routines in NASTRAN was emphasized during the 
orbiter analysis. New routines should be developed for NASTRAN which would generate 


31 



APPENDIX B 


all eigenvalues and up to about 50 eigenvectors for matrix orders up to approximately 
500 in reasonable times. Since new routines were not available, two of the existing 
eigenvalue extraction methods (see ref. 32, section 10) in NASTRAN Level 15.5 were 
used to obtain the natural frequencies and mode shapes for the combined orbiter in the 
substructuring procedure described previously and also for the individual substructures. 

Eigenvalue extraction by the Givens method . - The Givens method for extracting 
eigenvalues (transformation method based on the tridiagonalization techniques of Givens) 
was used for some of the phase I runs. Intermittent unscheduled exits occurred in some 
of these runs. Until errors are corrected, the Givens method in NASTRAN Level 15.5 
is considered unreliable for eigenvalue extraction. 

Eigenvalue extraction by the inverse power method with shifts .- The inverse power 
method with shifts for extracting eigenvalues was used on the phase II runs and some of 
the phase I runs. This tracking method is time consuming for large problems. Care 
must be exercised to ensure that absolutely all eigenvalues in the required range have 
been obtained. It is noted that a smaller field is required for this method in the READ 
(eigenvalue extraction) module than for execution of the SMPl module which performs 
the Guyan reduction. Also, less field is required for the inverse power method than for 
the Givens method for eigenvalue extraction. One or more eigenvalues and eigenvectors 
can be obtained on an initial or restart run which is checkpointed. As previously men- 
tioned, additional runs with modified frequency search ranges can be submitted until all 
the desired frequencies and modes are obtained. 

Restarts 

The restart capability in NASTRAN for utilizing data from previous executions was 
used on many of the orbiter runs. When using the appropriate Old Problem Tape (OPTP) 
on a subsequent run, the NASTRAN Data deck contained the Executive deck which included 
restart dictionary cards and an alter packet, a Case Control deck, and only modifications, 
if any, to the Bulk Data deck. The Bulk Data deck from the previously checkpointed run 
is contained on the OPTP. These cards, with two exceptions, must not be included again 
on a restart run since NASTRAN will abort the job if duplicate Bulk Data cards are found. 
However, the BEGIN BULK and ENDDATA cards - the two exceptions - must be Included. 

Unmodified restarts .- The NASTRAN Data deck for an unmodified restart included 
the full restart dictionary obtained on the previous checkpointed run and contained no 
modifications to the Bulk Data deck. The cards specifying DMAP alters to the rigid 
format prior to the last reentry point contained in the restart dictionary were removed 
from the alter packet in the Executive deck. 

Modified restarts .- Restart tables are included in each rigid format description in 
reference 25, section 3. These tables were useful in determining which modules had to 
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be reexecuted when modifications were made to the bulk data. Although many different 
modified restarts were made, only two of those most frequently used are described. 

When it was determined that the stiffness matrix had to be regenerated in the 
SMAl module, only the first card from the restart dictionary was used and the entire 
alter packet was included. This mode of operation was convenient for the users and 
ensured that all modules would be reexecuted. 

Substantial computer time was saved by executing modified restart runs when only 
the Real Eigenvalue Extraction Data card (EIGR) was changed. On these runs, the restart 
dictionary included the RESTART card with continuation cards through the last file check- 
pointed for reentry at the DPD module, which extracts the Eigenvalue Extraction Data 
from the Dynamics Data block. The DMAP sequence number for module DPD is 85 in 
Rigid Format 3. All DMAP alter cards specifying changes to the DMAP sequence prior 
to statement 85 were removed from the alter packet. 

Pseudo restarts .- Pseudo restart runs are runs to obtain only additional output data 
from calculations made on previous executions. Several restart runs were made durii^ 
the orbiter analysis. DMAP alters were inserted in the Executive deck to obtain print- 
outs for selected matrices. Modifications were made in the Case Control deck for print- 
ing and/or plotting as required. 

Changes in the plot package in the Case Control deck on a pseudo restart run will 
trigger execution or reexecution of only those modules related to plotting if a fully check- 
pointed run has been made previously. This mode of operation is desirable especially 
when solvir^ large problems. The primary reason for using this method to obtain 
deformed plots is that the field length required for plotting is substantially less than that 
required for solution of medium or large problems. Thus, several runs can be submitted 
with different plot packages for a nominal cost compared with the cost for a large solu- 
tion run. 


Plotting 

A post-processor system is utilized now for all graphics at the LRC computer 
complex. Use of this post-processor system, however, was optional during most of the 
orbiter runs. It was selected primarily to utilize on-line plotting capability which is both 
faster and less expensive than plotting on off-line devices available at the LRC computer 
complex. The general-purpose plotter was selected in the plot package in the Case 
Control deck. A PLT2 physical tape reel must be requested and mounted prior to 
NASTRAN execution whenever a plot package is included in the Case Control deck. 

In order to obtain NASTRAN plots at the LRC computer complex, three programs 
must be executed in order but not necessarily on the same job submittal. These programs 
are as follows: 
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(1) NASTRAN, which generates plot commands for the general purpose plotter on 
the PLT2 reel (It is noted that the PLT2 reel is not released by the LRC version of 
NASTRAN Level 15.5.), 

(2) NASTPLT, which reads the plot commands from the PLT2 reel and reformats 
them on the plot vector file, SAVPLT, in the format used by the LRC plotting post pro- 
cessors, and 

(3) PLOT, which formats the plot vector file from the SAVPLT file into a plot 
vector file for the particular graphic device selected on its load-execute card and also 
allows modifications to be made to the plot vector file. 

When on-line plots were obtained on the same job followii^ NASTRAN execution, 
the field length held by the job was reduced to 52 000 octal in order to release excess 
core to the operating system. 


Concluding Remarks 

NASTRAN is an effective and versatile tool for predicting vibration response of 
complicated structures. However, orbiter analyses at the Langley Research Center 
required long computer execution times with large field lengths, which, under normal 
operating procedures, resulted in extended turnaround times. These turnaround times 
could be shortened only by obtaining priority of operation in the LRC computer complex. 
This situation is equivalent to havii^ a dedicated computer available. The single DMAP 
alter packet and DMAP tape copy program provided by the Grumman Aerospace 
Corporation (refs. 21 and 26) for substructuring in NASTRAN is easier to use than the 
substructure analysis procedure described in the NASTRAN User's Manual (ref. 25). 
New eigensolution routines should be developed for NASTRAN which would generate all 
eigenvalues and up to about 50 eigenvectors for matrix orders up to approximately 500 
in reasonable times. 
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FLEXIBILITY INVESTIGATION OF ORBITER JOINTS 
AND DOOR LONGERON 

Arthur I. Miller* and Murray Bernstein* 

This appendix contains a brief description of the simulation of certain flexibilities 
not accounted for in the initial NASTRAN model. The most significant is the flexibility 
of the interface between the fin and the fuselage. Of the three models shown in figure 24, 
model FI was based on an approximate structural analysis of the fin-fuselage interface, 
whereas models F2 and F3 were derived from comparisons of static analytical and test 
deflection data. The use of model FI was limited to static analysis, but models F2 and F3 
were used in both static and dynamic analyses. Of lesser importance are the flexibilities 
of the juncture between the cabin and cargo bay and the additional bending flexibility of 
the door longeron due to its discrete -point attachment to the longeron along the cargo bay. 

Joint Flexibilities 

Modification of the NASTRAN finite -element model to account for joint flexibility 
was most readily accomplished by adding spring elements. Values for the spring con- 
stants were derived either from a simple idealization of the structural characteristics 
of the model joints or from correlation between static test deflection data and analysis. 

Fin-fuselage interface .- The interface shown in figures 23(a) and 23(b) is statically 
indeterminant and relatively complex. Loads to the fuselage are transmitted through the 
forward and aft spars by clips and by the fin-fuselage angle connection along the inter- 
face. Figure 23(b) shows the aft spar-clip connection. Loads through the center spar are 
transmitted by a gusset plate (fig. 23(a)) designed to transmit fin pitch into the main fuse- 
lage engine bulkhead. In attempting an analytical representation, it was assumed that 
only the forward and aft clips had to be considered and other elements were not signifi- 
cant. Such a model is designated as FI in figure 24. 


*Grumman Aerospace Corporation, 
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The actual specifications of the clips in model FI are shown in sketch (C1-). 
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The idealization of the clip as a pin-roller attachment is an attempt to determine the 
flexibility of the clip itself. The spring constant Kpj of the idealized clip is deter- 
mined from a consideration of bending energy alone. A representation of the fin-clip 
support is given in sketch (C2) : 



For the clip idealization shown therein, the strain energy in bending is given by: 
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The total deflection ^total bending is 
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The bending moment in span f j is given by 
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The spring constant Kpj of the idealized clip is given by 
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Model F2 shows a single spring to be added to the NASTRAN model to account for 
the difference between analyses and test results. Only a forward spring was required 
because, as noted in the main body of this paper, the aft attachment in the fuselage had 
a negligible effect on the frequencies and the aft fuselage support frame, to which the 
clip is attached, is quite flexible. 
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It is assumed that the measxired fore-and-aft deflection of the fin is composed of 
two components. One of these, 6a,s’ entirely to fin fore-and-aft bending and is 

found from aimlysis to be small. The other deflection component A is assumed to be 
caused by joint rotation. The spring constant is determined solely for joint rotation and 
is given by the following set of relations shown in sketch (C3): 



Taking moments about the pivot yields 
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Thus, the spring constant 
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The spring constant used in the dynamic analysis (Kj'2,dyn determined 

from the product of the spring constant used in the static analysis and the ratio of 
analytical-to-experimental deflections. 
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Model F3 in figure 24 shows the roll springs used in the antisymmetric loading. 
The roll springs were also determined from a correlation between experiment and 
analysis as shown by equations (Cll) and (C12) and in sketch (C4) (based on ref. 21, 
vol. n). 
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The spring constant Kp 3 is given by 
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Again, the spring -constant value for the dynamic analysis in fig. 24) was 

determined by multiplying Kj .3 by the ratio of analytical-to-experimental static 
deflections. 
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Comparison of finite -element models RSI and RS2 with static test (see table X) 
indicates that fin model FI allows excessive flexibility of the interface between fin and 
fuselage. Fin model F2 reduces the flexibility of the fin-fuselage Interface through 
elimination of the aft clip spring. This is possible since the fuselage support frame 
(which has no underlying frame structure) to which the aft clip is attached is already 
quite flexible, as is evident in figure 23(b). Other rod-shear panel models RS4 and RS5 
and reduced membrane model RMl yield fin deflections of sufficient accuracy (approxi- 
mately ±10 percent of measured fin deflection) when using fin model F2. Fin model F3, 
used in the case of antisymmetric loading, yields fin deflections within 5 percent of mea- 
sured data, as indicated by the results shown in table XI. 

The dynamic effects of joint flexibility are shown in table XII in a comparison of 
symmetric -mode orbiter frequencies calculated using model F2 flexibility with initial 
analytical frequencies and with measured frequencies. As may be seen, this modifica- 
tion resulted in about a 15- to 25 -percent discrepancy between analytical and measured 
frequencies. Fin-fuselage flexibility thus accounts for 13^ percent of the disagreement 
between initial analytical and measured frequencies in the first mode and for as much 
as 58 percent in the fourth mode. 

Cabin-to -payload bay juncture .- The CPB joint experiences primarily axial loadir^ 
during fuselage vertical bending. The axial loading is assumed to be carried by the 
aluminum splices across the 0.0813-cm (0.032-in.) joint gap. The upper splice is fur- 
ther assumed to transmit bending due to the 6° bend. The upper splice may be idealized 
as a pin-roller configuration for bending deflection as in model FI of figure 24. 

Joint flexibility for the upper splice is based on the deflection relation 

^upper = ^x,axial + ^x, bending (C14) 

and for the lower splice 

®lower " ^x, axial (C15) 

The structural idealization of the upper splice is illustrated in sketch (C5): 
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Then, for the upper splice, 
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and the spring for the lower splice is given by Vj ^lower’ 

1 ApE 

" bp/ApE bp 


(C21) 


Assuming that the upper and lower splices act as springs in parallel, a new sprii^ con- 
stant for the complete joint is simply the sum of the two spring constants K^ and K^. 
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The value of the total spring constant is given in figure 25. The greater portion of 
this value comes from since no bending deflection is included as it is in the value 
for Ky. The actual values of and Ky are given as follows: 

= 28 895 kN/m (165 000 Ib/in.) 

°x, axial 

Ky = = 1999 kN/m (11 412 Ib/in.) 

^Xjaxial °x, bending 

Assessment of the CPB spring on static deflection cannot be made since no attempt 
was made to isolate its effect from other aspects of the analytical models. Dynamically, 
the CPB spring had negligible effect on frequencies and mode shapes in fuselage phase I 
(substructure) analysis. 


Door Longeron Effectiveness 

In the initial NASTRAN model, the door longeron was considered fully effective 
during fuselage bending. However, the door is attached physically to the fuselage by 
four shear pins spaced approximately 63.50 cm (25 in.) apart. One pin is shown in the 
door clip of figure 6. These are fastened into clips which brii^ the restraint point down 
about 2.54 cm (1 in.) below the neutral axis of the door longeron. Any axial (x-direction) 
load applied by the shear pins will result in a combination of extension and rotation of the 
door longeron. The bending flexibility due to rotation had not been accounted for and was 
therefore added. The rotational flexibility of the longeron as a beam was calculated as 
a small separate NASTRAN problem and translated into an effective reduction in cross- 
sectional area from 0.3613 cm^ (o.056 in^) for the original longeron to 0.1290 cm^ 

(o.02 in^). This modification resulted in a decrease in calculated fuselage bending stiff- 
ness of about 4^ percent. An alternative approach adopted in reference 21 was remodel- 
ing this area by revising the geometry and adding additional grid points and constraint 
relationships. 

The calculated flexibilities (i.e., deflections for unit loads) with the doors on were 
a reasonable match with measured flexibilities, but with doors off, the analytical model 
was too flexible. The results of the static test, upon initial review, indicated no signifi- 
cant difference in deflection with and without the doors. However, strain gages on the 
door indicated that the door-longeron stresses were about half those in the fuselage 
longeron and therefore the door was partially effective. Dynamically, the reduction in 
longeron area caused a negligible effect on frequencies and mode shapes. 
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Concluding Remarks 

Investigation of the flexibility at the fin-fuselage interface has shown that fin 
models F2 and F3 yield fin deflections within 10 percent of measured deflections. 
Dynamically, fin model F2 accounts for approximately 13^ percent of the first-mode 
frequency discrepancy between initial analysis and test. The effect of the CPB sprii^ 
and reduction of the door loi^eron area result in negligible effects on overall orbiter 
frequencies and mode shapes. 
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NASTRAN MODELING ACCURACY STUDIES 
Jerrold M. Housner and Manuel Stein 

The purpose of this appendix is to present and discuss results of^two studies which 
were undertaken to evaluate NASTRAN modelii^ accuracy for dynamic response of the 
fuselage of the 1/8-scale shuttle orbiter model. In the first study, the general capability 
of NASTRAN to predict accurately the response of the fuselage -like structure shown in fig 
ure 41 is considered, where shear lag and cross-sectional warping are included. The sec 
ond study is concerned with the effect on the natural frequencies and mode shapes of high 
aspect ratio of NASTRAN’s QDMEM2 membrane elements used in the orbiter fuselage. 

Symbols 

aj,a 2 jag amplitudes of deflection functions for symmetric modes 
aj,a. 2 ,a 2 ,a^ amplitudes of deflection functions for antisymmetric modes 
Aj,A 2 ,A 3 ,A 4 cross-sectional areas of fuselage stiffening rods 
AjjA^jAjj cross-sectional areas of floor, wall, and keel, respectively 
Afj. total cross-sectional area, 2^Aj^ + A 2 ) + Ag + A^ + A^ + 2A^ + Aj^ 

b semichord of fuselage cross section 

Q 

9 warping coefficient given by equation (D18) 

Cj coefficient given by equation (D7) 

E Young's modulus 

fy coefficients defined by equation (Dll) 

G shear modulus 

h-y^,hjj heights of wall and keel, respectively 
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moments of inertia about axes passing through the cross-sectional 
centroid and parallel to the Y- and Z-axes, respectively 

polar moment of inertia of cross section about the origin 

effective moment of inertia defined by equation (D5) 

torsion constant given by equation (D17) 

length of simplified fuselage 

mass per unit length of the fuselage 

cross-sectional coordinate, as shown in figure 42 

thicknesses of keel, floor, and wall, respectively 

axial displacement (x-direction) 

transverse displacement (y-direction) 

lateral displacement (z -direction) 

deflection functions defined by equations (Dl) to (D3) 

Cartesian coordinate system (see fig. 42) 

shear -center offset from centroid in z-direction 

shear -center offset from Y-axis 

clockwise rotation of cross section, as well as displacement function 
defined by equation (D3) 

displacement fxmction defined by equation (Dl) 
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V Poisson's ratio 

u) angular frequency, rad/sec 

p mass density 

The primes indicate differentiation with respect to x. 

Methods of Analysis 

To confirm the NASTRAN results for a simplified fuselage model, two closed-form 
mathematical analyses were performed. One analysis, developed herein, was based on 
an energy approach usirg an assumed displacement state, piecewise linear along the 
cross-sectional perimeter. The second analysis was based on a thin-walled, open- 
section theory (ref. 33). Included also are descriptions of two NASTRAN models of the 
simplified fuselage of figure 41. In contrast to the l/8-scale orbiter fuselage, the sim- 
plified fuselage has uniform geometric and material properties along its length. 

Piecewise linear analysis .- Due to symmetry, only half the fuselage cross section 
needs to be considered, as shown in figure 42. This model may be viewed as consisting 
of three plates - the keel, the floor, and the wall - and periodically spaced frames. The 
fuselage frames (fig. 41) are considered to be stiff enough so that the strain €g and cur- 
vature Kgs in each plate is zero. The subscripts here are associated with the cross- 
sectional coordinates running counterclockwise around the cross section as shown in fig- 
ure 42. Under this stiff -frame assumption, the deformations of the middle plane of each 
plate may be expressed as the following piecewise linear functions of s: 


Uj(x,s) = Uj(x) + Sp^ix) 

(Dl) 

Vj(x,s) = Vj(x) 

(D2) 

Wj(x,s) = Wj(x) + S0j(s) 

(D3) 


where i = 1, 2, and 3 identifies keel, floor, and wall, respectively. In equation (Dl) the 
inclusion of the function x(/i{x) allows for the warping of the cross section, while in equa- 
tion (D3) the function 0^(x) allows for the rotation of the cross section. 

Equations (Dl) to (D3) indicate 15 unknown functions, but not all of these are inde- 
pendent. Since the displacements must be continuous at the joints and right angles must 
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be preserved at the joints, and dw^^ds must be continuous. This implies eight 
continuity constraints on the unknown functions. 

Symmetric modes: For symmetric modes, the following additional constraints may 
be imposed: 




(D4) 


Consequently, upon application of the continuity and symmetry constraints, only three 
unknowns remain in equations (Dl) to (D3). Calculation of the strains in each plate and 
application of the variational principle finally leads to the following equations of motion: 


\ 2 2A^G ^ 


o_ 

mco^Wr 




eI 


(D5) 


(D6) 


u 


2 


II 




ft 


where 


(D7) 


I = (hjA^ + hk^Ak) + 2Aih^2 ^ 


3(1 - v^) 


-Cl 


1 _ ^ 2 !^^^ ^ ^k^k) -^4^k| 


and 


^wAw + hkAk + (2h^Aj + A4hk)(l - 

Cj ■ r 

Af + 2A^ + Ak + (2Aj + 2A2 + A3 + A4)(1 - v^) 

The solution to equations (D5) to (D7) for simply supported boundary conditions is 


W 2 - sin 

(D8a) 

= ^2 cos 

(D8b) 

Ug = ag cos 

(D8c) 
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Substitution of equation (D8a) into equation (D5) yields the frequency relation 

n^ 

El/m 1 , Ef (m\^ 


When — S' « 1, the lowest frequency is found for n = 1. Once the frequency is 

2GA^je^ 

determined, the corresponding mode shapes may be obtained from equations (D5) to (D8). 


2GA^je^ 


Antisymmetric modes: For antisymmetric modes, the following constraints 
replace equation (D4): 

ui = il/i = vi = 0 (DIO) 

Consequently, upon application of continuity and antisymmetric constraints, there are 
four unknowns remaining in equations (Dl) to (D3). Calculation of the strains in each 
plate and application of the variational principle result in the following set of equations: 


(Dll) 


where the elements of matrix []fj are 


- p:.2[2b2A„ + I + 1 i 


+ 2Ij^ + 41^ + 




+ 2b2(Aj + Ag) + 2hJ^A^ + hjj.^ + Ij] 

f 12 = ^21 = (xf ^ I Vk) - I l^k\ ^ 2h^^l 


^13 “^31 “ ° 


f 14 = f^i = -2bGiV ^ 
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Nontrivial solutions to equation (Dll) are obtained from the determinantal fre- 
quency equation 


1 f 1= 0 

(D12) 

which yields eigenvalues and eigenvectors. From the eigenvectors 
may be found from the following relations: 

{a^, the mode shapes 

01 = 02 = ^3 = ai sin ^ 

(D13a) 

?2 = ag sin 

(D13b) 

^2 = -h T 

(D13c) 

;^^3 = a4COsSf 

(D13d) 


Analysis of thin-walled open section .- In reference 33, Chajes and Winter present 
a theory for the buckling of thin-walled open-section columns. In order to make this 
theory applicable to vibrations, the buckling terms in their equations are replaced by 
the appropriate inertial terms; and advantage is taken of symmetry to reduce the equa- 
tions to 
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ElyW^ - mw^w = 0 


ElgV^ - mco^v - mco^ZQ0 = 0 


TV T t Iq 

EC00^ - GJ0 




9 - mio^ZQV = 0 


(D14) 

(D15) 

(D16) 


where 


J = + V 3 ^) 


Ce = I b^t2[e3 - (h„ - £f] + I eV® + 2Aib2(h^ - ef + 


(D17) 

(D18) 


The term ECg0^ provides the inherent resistance of the cross section to warping; 
ly and I 2 are the cross-sectional moments of inertia about axes passing through the 
centroid and parallel to the Y- and Z-axes, respectively; Iq is the polar moment of 
inertia of the cross section about the origin; Zq is the offset of the shear center from 
the centroid in the z -direction; e is the offset of the shear center from the Y-axis. 

Symmetric modes: For modes which are symmetric about the plane y = 0, 

0 = 0 and vibrations are governed by equation (D14), which is clearly the simple beam 
equation without transverse -shear effects. 

Antisymmetric modes: For modes which are antisymmetric about the plane y = 0, 
w s 0 and V and 9 are given by 


V = V sin 


nirx 


(D19) 


0 = 0 sin 


nTTX 


(D20) 


Substitution of equations (D19) and (D20) into equations (D15) and (D16) yields the follow- 
ii^ frequency equation for nontrivial solutions: 


El, 


4 ¥) - 


-mctj^z^ 


9 

-mco z. 




= 0 


(D21) 


The lowest frequency of equation (D21) occurs for n = 1. 
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NASTRAN models .- Two NASTRAN models of the simplified fuselage of figure 41 
were considered and are shown in figure 43. As a consequence of symmetry, only a 
qviarter of the fuselage needs to be retained. Both coarse (fig. 43(a)) and refined 
(fig. 43(b)) models were necessary to confirm convergence of solutions. Before Guyan 
reduction the coarse and refined models for the quarter fuselage contained 128 and 
746 degrees of freedom, respectively, for the symmetric modes, and 112 and 694 degrees 
of freedom, respectively, for the antisymmetric modes. Comparison of these numbers 
(doubled) with those in table V shows that the coarse and refined models bracket the 
1/8- scale fuselage model, with the refined model having nearly the same number of 
degrees of freedom as the 1/8 -scale fuselage model. 

A simple support at the end x = 0 was provided by setting v = w = 0 at all 
nodal (or grid) stations around the end. Setting the longitudinal displacement u, 
together with pitch and yaw rotations, equal to zero at all grid stations at the fuselage 
midspan gave the symmetric modes which included the lowest modes of interest. 

For modes which are symmetric about the plane y = 0, the longitudinal displace- 
ment u, together with roll and yaw rotations, were set equal to zero at all grid points 
in the plane. For modes which are antisymmetric about this plane, the longitudinal 
displacement u, the vertical displacement w, and the pitch rotation were set equal to 
zero at all grid points in the plane. 

Evaluation of NASTRAN Models for Vibration Analysis 

The general capability of NASTRAN to predict accurately frequencies and mode 
shapes of the fuselage was determined by compariig NASTRAN solutions with the closed- 
form mathematical solutions previously described for the simplified model of figure 41. 
The geometric and material properties were considered uniform along the length of the 
model and are listed in table Xfll. The length of the model was chosen to correspond 
to the distance between nodal points of the 1/8-scale fuselage in its first symmetric or 
antisymmetric mode, and, as previously noted, the ends were taken to be simply 
supported. As also previously observed, fuselage symmetry about the plane y = 0 
enables the separation of modes into symmetric or antisymmetric about this plane. 
Furthermore, for the antisymmetric modes, the shear center - though lyirg in the 
plane y = 0 — does not coincide with the cross-sectional centroid. Consequently, 
lateral bending and torsion motions are coupled. 

Frequencies calculated by both closed-form mathematical analyses are compared 
with NASTRAN frequencies in table XTV for the first symmetric and antisymmetric modes. 
In general, the results for each mode are in good agreement, and in particular, the refined 
NASTRAN model (fig. 43(b)) is in excellent agreement with the piecewise linear analysis. 
For the symmetric mode, the frequency by the beam solution is 12 percent higher than 
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that by the piecewise linear analysis. As is evident from equation (D9), the sole dif- 
ference between the two solutions is the inclusion of transverse shear in the piecewise 
linear analysis. Hence, transverse shear is important in this mode. 

In figure 44, the axial displacements and shearing stresses at a simple support 
are shown for the NASTRAN refined model and the piecewise linear model for the first 
symmetric mode. The agreement is seen to be excellent. In figure 45, the axial dis- 
placements and shearing stresses at a simple support, along with the lateral and trans- 
verse displacements at the fuselage midspan, are shown for the first antisymmetric mode. 
Again the agreement in displacements is excellent. The agreement in shearing stresses 
is also very good; however, as is often the case, the agreement is not as good as that for 
the displacements. It is also of interest to note in figure 45(b) the sizable amount of 
cross-sectional warping which takes place. 

Although not shown, the modal displacements predicted by the coarse model were 
up to 6 percent less than those of the refined model. It is not reasonable to compare 
the modal stresses predicted by the coarse and refined models, since the NASTRAN pre- 
dicted stresses are not at specific points, but are averages over an element, and the 
coarse elements are about six times bigger than the refined elements. 

In summary, the coarse NASTRAN model (fig. 43(a)) was within 6 percent of the 

piecewise-linear solution, whereas the refined model provided excellent results (within 

1 \ ^ 

1^ percent of the piecewise-linear solution). 

Aspect Ratio Characteristics of NASTRAN Membrane Elements 

The NASTRAN model of the 1/8-scale orbiter fuselage is largely composed of flat 
quadrilateral membrane finite elements. In general, these elements are too stiff in 
in-plane bendii^ and must therefore be used cautiously in portions of a structure sub- 
jected to in-plane bending (see ref. 34). In particular, the QDMEM2 membrane element, 
used in the NASTRAN model, loses significant accuracy at high aspect ratios. This was 
demonstrated in an unreported analysis by William C. Walton, Jr., and Huey D. Carden 
of NASA Langley Research Center. Using a cantilever beam made up predominantly of 
QDMEM2 elements of aspect ratio 5 and subjected to a static concentrated load at the 
beam tip, they obtained an in-plane bending stiffness that was about 37 percent higher 
than the classical beam solution. This increase corresponds to a frequency 17 percent 
too high. With panels of aspect ratio 3, the stiffness increase was 17 percent, which 
corresponds to an 8-percent frequency increase. 

Since the QDMEM2 elements were used with aspect ratios as high as 4 in the 
1/8-scale orbiter model, it was considered necessary to examine their effect on NASTRAN 
predicted fuselage frequencies. Three separate studies were undertaken. 
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In the first study (ref. 21, vol. II) by Philip W. Mason (GAC), it was found that in 
modeling a structure like the fuselage with spar caps on the tension and compression 
sides, the stiffness increase due to high-aspect-ratio QDMEM2 elements is reduced as 
the spar-cap cross-sectional area is increased. Relating this result to a typical fuse- 
lage station. Mason estimated the frequency increase would be less than 5 percent. It 
is reasonable to conclude that when the spars are stiff enough to carry the burden of the 
bending loads, the membrane elements, acting together, behave like a shear web, and, 
since the membrane elements are designed to resist in-plane shear, the error due to 
their presence is reduced. 

In an unreported second study by John L. Sewall (LRC), the vertical membrane 
fuselage sidewall in the payload-bay area, forward of the wing carry -through structure, 
was remodeled with an equivalent rod- shear panel combination. Imposing in-plane 
static loading on both the membrane and equivalent sidewall models resulted in only a 
2- to 3-percent increase in stiffness of the membrane panel over that of the equivalent 
rod-shear panel model. Further examination showed that this difference for the side- 
wall could be attributed to interspersing panels of aspect ratio 1 and 2 among panels of 
aspect ratio 3 and 4 in the original membrane model. 

Finally, in the third study (reported herein), a NASTRAN model of the simplified 
fuselage shown in figure 46 was considered. This model was composed of nearly all 
aspect-ratio-3 QDMEM2 elements; and an attempt was made to have, on the average, 
about the same number of elements in both the sidewalls and the floor of the fuselage 
as was in the original 1/8-scale NASTRAN fuselage model. This study showed only 
a 2-percent frequency increase compared with the refined NASTRAN model of fig- 
ure 43(b), in which the panel elements were of aspect ratio 1. 

In summary, it appears that the use of relatively high aspect-ratio QDMEM2 ele- 
ments in the orbiter NASTRAN model will raise the frequency only a few percent. This 
increase would probably be higher if all the elements were of high aspect ratio or if 
spars were not present in the fuselage. 


Conclusions 

Results of the investigation reported warrant the following conclusions: 

1. In general, NASTRAN solutions of a simplified orbiter fuselage structure for 
vibration mode shapes and frequencies correlated reasonably well with solutions by 
closed-form mathematical analyses. A coarse NASTRAN model with two panel ele- 
ments in the walls and four in the floor gave frequencies within 6 percent. Refining 
the model to Include 5 elements in the walls and 10 in the base greatly improved the 
accuracy to within 1^ percent. 
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2. The use of relatively high aspect-ratio membrane elements (QDMEM2) in the 
or biter NASTRAN model does not result in a frequency increase of more than a few per- 
cent. This increase is expected to be higher with aD panel elements of high aspect ratio 
and/or without spars on the tension and compression sides of the fuselage. 
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PANEL IMPERFECTION STUDIES 
Murray Bernstein* and John L. Sewall, Coordinators 

In the original NASTRAN model, membrane elements makii^ up the external fuse- 
lage and wing surfaces were considered fully effective in transmitting in-plane direct and 
shear stresses. However, inspection of the physical model showed that most surface 
panels were not perfectly flat but had initial out-of-plane bow considerably more in mag- 
nitude than the panel thickness. This deviation suggested that the panels could not be 
fully effective in carryii^ the in-plane stress. This was especially true of panels having 
large unsupported areas. This ineffectiveness was further evident in observations of 
localized oscillations durii^ vibration testing and in strain-gage readings from static 
tests. To obtain an indication of the sensitivity of the vibration modes to reduced panel 
effectiveness, the vertical sidewalls of the analytical fuselage were remodeled with rods 
and shear panels replacing the membrane panels, and a complete NASTRAN vibration 
analysis (i.e., through phase II, see fig. 21) was made. The rods were of token or 
minimal, cross-sectional areas to stabilize the shear panels. The results of this analy- 
sis appear in table XV. 

The significant frequency reductions produced by this modification stimulated a 
search for a valid measure of panel effectiveness. This search was followed along two 
general paths: First, static test deflections were compared with static analytical deflec- 
tions to determine and evaluate candidate panel modifications; and second, independent 
studies were conducted both at LRC and GAG to determine effectiveness criteria. 

Results from both approaches were applied in NASTRAN vibration analyses for 
comparison with measured frequencies and mode shapes. Frequencies of all the panel- 
effectiveness models are listed in table XVI aloi^ with the initial analytical frequencies 
and measured frequencies. 

PANEL EFFECTIVENESS BASED ON STATIC TEST DATA 
Arthur I. Miller*, Jack R. Barrett**, and Murray Bernstein* 

Panel Effectiveness in Direct Stress 

The initial analytical deflection ctirves when compared to static test points showed 
a uniformly stiffer appearance than the test model deflection curve (see figs. 47 to 49). 

Grumman Aerospace Corporation. 

** Rockwell International Corporation. 
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The continuity of the deflection curve suggested a general overstiffness rather than a 
localized overstiffness. A general overstiffness was suspected in modeling the panels 
with fully effective skins. In modeling beam webs by rods and shear panels, one-sixth 
of the panel cross-sectional area is concentrated in the cap areas to represent a fully 
effective web in bending (as shown in ref. 35). However, the panel shear capability 
remains fully effective. This concept is applicable in representing thb fuselage side- 
wall in vertical bending and the fuselage bottom panels in lateral bending. 

When modeling skin-stringer structures subjected to compressive loads, the 
stringer and some skin is assumed effective. This is the case with wing cover sheets 
in wing bending and fuselage sidewalls and bottom panels in lateral and vertical bend- 
ing, respectively. The problem then is to determine how much skin is actually effective. 
A number of analyses with varying amounts of effective skin could be run in order to 
achieve a correlation between static test deflections and analysis. However, it was 
considered more expedient to solve for an amount of effective skin in the test model 
directly and use it in the revised analytical model. 

From the deflection equation normalized to a unit load 

5 oc F (F = 4.45 N (1 lb)) (El) 

where the constant of proportionality involves length cubed divided by bending stiff- 
ness EL The modulus E and lengths of the analytical model were considered to 
be equivalent to those of the test model, thus allowing only the flexural moment of 
inertia I to vary by the amount of effective skin. The equation relating test and 
analysis can be stated as 

6tIt=6aIa (E2) 

and the moment of inertia of the test model I-p can be expressed as the ratio of maxi- 
mum analytical-to-test deflections niultiplied by the moment of inertia for 

the initial analysis I^, based on fully effective panels; i.e.. 



This effective I-p in the test model now becomes the required moment of inertia in the 
revised analytical model. The skin area required in the analytical model to satisfy I-j* 
is formulated for the fuselage and wing in the next two sections. 

Fuselage representation.- The fuselage is idealized as shown in sketch (El). 
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, denote existing 
longeron areas 


Fuselage effective skin areas A-j., Ag, and A^ are determined by assuming that the 
neutral ^is location of a fuselage section with fully effective skins is maintained, even 
though the moment of inertia for the section is reduced. The neutral axis of a cross 
section in terms of the unknown areas A-p, Ag, and Aq is given by 


z = 


1 

jAj + A^j 

D + A^C 

(Ai + A2 + A3 + A^^ 

1+ 1 

[At + Ag + Aq) 


(E4) 


Moments of inertia required at a cross section are determined by the following: 



T/Z 


^A,Y - (^1 + A-j. 


vr2 


i^D - z)^ + ^A2 + A3 + Ag + Aq^z^ + A4(C - z)' 


(E5) 


and 



where 


(E6) 



ratio of maximum deflections between analysis and test in the z -direction 


ratio of maximum deflections between analysis and test in the y -direction 


and the moments of inertia of the rods about their own centroids are neglected. 


The terms 2 ^A Y ^®P^®sent moments of inertia of fuselage cross 

sections with fully effective skins about the neutral Z- and Y-axes. Solution of equa- 
tions (E4) to (E6) yields the three skin areas. These areas are: 
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Afji = 


ly - AiD(D - z) + A 4 C(C - z) 
D(D-z) 


(E7) 


A 


B 


Iy-(Ai + At)(D-^2_a4(C-z)2 ^ ^ ^ 

Zn (A2 + A3 + Aq) 


(E8) 


Ac 


ly + A 4 C(C + D - 2^ 
z(D - z) 



(E9) 


The values of these areas over representative fuselage stations, as shown in table XVII, 
indicated that they were equivalent to having 46 percent of the cross-sectional skin area 
in the cargo bay effective in direct stress. The three areas were each represented by 
rod elements at the indicated locations adjacent to existing longerons, and the skin was 
modeled by shear panels, together with minimal rods where actual rods did not exist. 

The NASTRAN model RS2 was modeled in this fashion. 

Included also in table XVII are effective areas arbitrarily distributed according 
to the 1/6 -area criterion described earlier and amounting to 33 percent of the total skin 
cross-sectional area. This criterion was used to model RSI. 

Wing representation . - The effective skin for carrying direct stress in the wing and 
wing carry-through structure was found in a manner similar to that for the fuselage. The 
total area moment of inertia Ia,X initial analytical model at various cross sec- 

tions alor^ the span was reduced by the wing-tip deflection ratio 
required inertias Ix for the revised analytical model as follows: 


ix = 



(ElO) 


The deflections involved in this equation are shown in the following sketch (E2): 



Sketch (E2) 
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A simplified typical wing cross section between two spars is represented in sketch (E3). 

Upper 
Lower 

Sketch (E3) 

The required moment of inertia Ix for the revised analytical model consists of two 
parts 

Ix = Ic + Is (Ell) 

where 1^, is the moment of inertia associated with in-plane bending of the spars about 
the neutral axis 

Ic = AiZu^ -I- 



and where Ig is the moment of inertia of the wing cover skins effective in compression 

Is = w(t^Zu^ + h'^^) 


The areas Aj and A 2 contain the top and bottom spar -cap areas plus their associated 
1/6 web areas (in accordance with ref. 35), and W is the effective skin width ^ ^i 
where Wj is the effective skin width to one side of a spar. With Ix obtained from 
equation (ElO) and I^, known from cross-sectional geometry, the only unknown equa- 
tion (Ell) is W, which may be expressed as 


Iy - Ir 




(E12) 


Knowing W and the number of caps present in the cross section, Wj for each cap can 
be found and the reduced skin area Wj^t can be input into the analysis. The portion of 
skin effective in a section is given by 


W _ ^uW t^W A 

Wp tu^F ^total 


(E13) 


where A is the effective skin area and is the total skin area. 
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For the spanwise stations checked, the percent of skin effective was found to vary 
only slightly from 76.5 to 81.3 and averaged 80 percent. The wing was then modeled 
using shear panels and rods which included 80 percent of the skin areas and the total 
spar -cap areas. Minimal rods were again used for stability around panel e<%es not 
bounded by existing rib- or spar -cap rods. Rib-cap areas were arbitrarily sized using 
the 80-percent ratio. This method was used on model RS2 and gave the best static deflec- 
tion correlation with test, as shown in figure 49. 

The panel effectiveness in direct stress for the wing is indicated for models RSI 
and RS2 in table X in the wing -tip row. In model RSI, only the carry -through struc- 
ture panels were reduced for the wing, and the panels outboard of the fuselage were 
fully effective as in the original analysis. As in the case of the revised fuselage for this 
model, the wing carry-through panels were also arbitrarily modeled by the 1/6 -area 
criterion of reference 35, with 1/6 of the web areas lumped in the top-and-bottom spar- 
cap areas for a 33 -percent panel effectiveness. The deflection curve for this model in 
figure 49 shows that the desired wirg-tip deflection was obtained with only this inboard 
part of the wing structure reduced in stiffness. Since both outer wing and its carry - 
through structure were of the same type of construction, a uniform reduction was needed 
for both parts to be more effective than 33 percent, i.e., if 1/6 of the panel area had been 
used for the outer wing as well as for the carry through, the desired deflection would 
have exceeded the test value. With available static test data, it was possible to find the 
percent skin effective in the test model and verify this percentage in an analytical model. 
This was done for model RS2, as described in the foregoing paragraphs, and resulted in 
a tip deflection within 2 percent of the test value. 

Model Evaluation in Direct Stress 

Table X provides a summary of static deflection results for several analytical sym- 
metric models. The degree of flexibility of these models over the original is apparent 
from the table, but, it is difficult to judge from these data alone as to which of the two 
models, RSI or RS2, yields the best deflection results. However, the static deflection 
curves in figures 47 to 49 show that correlation between static test and analysis of the 
fuselage and wing was best for analytical model RS2. 

Panel Effectiveness in Shear 

In addition to investigatiig the effectiveness of flat imperfect panels in direct stress, 
modifications were also introduced in an effort to accoimt for reduced in-plane shear effec- 
tiveness. This was accomplished in three steps. The first two steps apply to rod- shear 
panel models whereas the third step involves a reduced membrane panel model. 
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In the first step, it was assumed that the effectiveness in shear would be the same 
as in direct stress. The resulting effect on static deformation is noted by comparing 
models RS4 with RS2 in table X. All panels of RS2 are fully effective in shear, whereas 
in model RS4 panel thicknesses are reduced so that the shear effectiveness was 46 per- 
cent for the fuselage cargo -bay region and 80 percent for the wii^. Deflections for loads 
at the midfuselage increased 20 percent due to added shear effects; however, the increase 
was less than 10 percent for other loading conditions. The fin deflection should not be 
considered in the comparison because of the chaises in the fin-root flexibility represen- 
tation from FI to F2. Vibration analysis usii^. model RS4 resulted in eigenvalues much 
closer to the measured values than the original model but still somewhat high (see 
table XVI(a). 

As a second step in evaluating effectiveness in shear, all panels includii^ those in 
the fuselage were assumed 80-percent effective. The resulting deflections are listed 
under model RS5 in table X. The midfuselage value is now within 5 percent of the corre- 
sponding value for model RS2. Similar comparisons for loads at the fuselage nose and 
cabin ballast are not valid because a required constraint relationship at the junction of 
the cabin and fuselage longeron which had been inadvertently omitted in previous models 
was now added. The resulting change in deflection was not considered large enough to 
require rerunning the large computer programs for the other models. The deflections 
along the length of the fuselage are shown in figure 48 for midfuselage load. Model RS5 
is slightly more flexible than the physical model (i.e., analytical deflections higher than 
measured deflections), and the shear effect is not significant for this static loading. 
Vibration analysis of this model resulted in frequencies still about 7 to 10 percent higher 
than measured values as shown in table XVI(a). 

In order to determine the separate sensitivities of the fuselage and wings to panel 
imperfections, a vibration analysis was made with fuselage skin in the payload-bay area 
fully effective in direct stress and shear but with the wing panels 80-percent effective. 
This model is designated RS3 in table XVI(a), and, as can be seen, there are negligible 
changes in all but the second frequency, which is a predominant wing-bendir^ frequency. 
For this mode, modifications to the wii^ panel effectiveness resulted in a 5-percent fre- 
quency reduction. Comparison of these frequencies with those of model RS5 indicates 
the wii^ to be far less affected by panel imperfections than the fuselage. 

The third step in accounting for in-plane shear effectiveness involved returning to 
the membrane panel representation and modifying it in accordance with results from 
separate panel effectiveness studies. In contrast to the rod-shear panel model, in-plane 
shear effectiveness in the membrane model cannot be different from direct-stress 
effectiveness. Work at GAC (ref. 21, vol. II) based on a nonlinear analysis indicated 
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that for panels typical of the fuselage with imperfections on the order of 1.75 times the 
panel thickness a direct stress effectiveness of about 0.55 would be representative of the 
model. Tests on a simplified panel model described in the subsequent section indicated 
an effectiveness of about 0.5. Additional studies at LRC based on an elastic finite -element 
analysis gave a value of 0.63 for panel effectiveness. It was decided that a direct stress 
effectiveness value of 0.6 would be representative of the panel aspect ratios and initial 
imperfections found in most places on the model. Moreover, as discussed in refer- 
ence 21, vol. II, the available information on shear indicated effectiveness values, both 
of 0.6 and 2/3. 

Model Evaluation in Direct Stress and Shear 

The 0.6 -effectiveness value was used for models RMl and RM2. Thicknesses of all 
exterior membrane panels in the fuselage, wings, and wing cover panels were uniformly 
reduced to 0.6 their original values. In model RM2, this thickness reduction was extended 
to membrane panels in the cargo-bay doors. The deflections of model RMl under most 
static loads were higher than measured values, as may be noted from figures 48 and 49 
and from table X. However, the model deflects less than the measured values at key 
significant locations; namely, the payload and cabin ballast where most of the weight is 
located. For that reason, the frequencies obtained from the vibration analysis shown in 
table VI are still higher than those measured. This difficulty of having the model too 
flexible for some static loads, yet higher infrequency, indicated that additional correc- 
tions were required. This observation can be drawn from comparison of the analytical 
and measured mode shapes in figures 12 to 15 and 26 to 29. However, since the revised 
model was considered close enough to the measured values to be adequate for appropri- 
ately combining the or biter with the external tank and solid rocket boosters (e.g., by 
some modal coupling) and because of other higher priority work involving these other 
two shuttle elements, it was decided to forego any further modifications. An independent 
effort at modifying the model was continued at GAC (see ref. 21) where some of the 
anomalies were resolved. 

The same membrane modifications were also made for the antisymmetric model, 
and static deflections for model RMl are compared with those for the original analytical 
model and with measured deflections in table XI for four antisymmetric loadings. Fre- 
quencies are compared with initial analytical and measured frequencies in table XVI(b) 
and are generally in better agreement with measured frequencies, as are those of 
model RS6. Extending the 60-percent panel effectiveness for model RMl to membrane 
panels in the door (model RM2) resulted in a negligible increase in flexibility and a 
negligible decrease in frequency. In general, the 0.6 -effectiveness factor introduced 
sufficient flexibility into the ordinal analytical model to significantly improve its agree- 
ment with the physical model but, in contrast to the symmetric model, not so much as 
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to make the analytical model more flexible than the physical model. This was true with 
doors on or off except for the midfuselage lateral load, for which the doors-on case of 
the revised analytical model was more flexible than the corresponding physical model. 

PANEL EFFECTIVENESS DETERMINED FROM INDEPENDENT STUDIES 
Murray Bernstein* and John L. Sewall, Coordinators 

Independent studies that guided the choice of the panel effectiveness factor for the 
reduced membrane models (RMl and RM2) are described in the following sections. These 
studies consisted of an analytical invest^ation at GAC and tests and analysis of a simple 
reinforced panel model at LRC. 


Theoretical Criteria 

The work reported in reference 21 involved calculating the effective panel width 
under direct stress loading for panels of aspect ratios and load intensities representa- 
tive of the fuselage and wing. Timoshenko's equations were used for plates simply 
supported along the side and axially loaded along the ends. The displacement functions 
and assumed initial bow were similar to those in reference 30, and the range of maximum 
displacement of the initial bow is extended to twice the panel thickness. The method of 
solution involved writing the expression for potential energy in terms of the displacements 
and initial imperfections, determining the minimum value of potential energy with respect 
to the displacements, and solvii^ for the strain and stress distributions. The stress dis- 
tribution was then integrated to get the panel load and effective width. Typical analytical 
results for three panels are shown in figure 50. 

These calculations showed that for panels with aspect ratios of 1.32, which are 
typical of the bottom of the fuselage, the effectiveness varies from 0.9 for maximum 
displacement of the initial bow of 1/4 thickness to 0.5 for twice the panel thickness at 
low load levels, as indicated from the inset plot of figure 50(a). Similar calculations 
for panels of 1.79 aspect ratio typical of the side wall, showed effectiveness to vary 
from 0.72 for 1/4 thickness initial bow to 0.47 for twice the thickness. The wing panels, 
which have aspect ratios of 0.54, were very effective even with maximum displacement 
of the initial bow up to twice the thickness. This work was not extended to shear and 
bending. 

However, work in references 36 and 37 indicated that an effectiveness of 2/3 in 
shear would be applicable. This was reinforced by data from previously unreported 
tests of riveted panels at GAC in which the data, although widely scattered, indicated an 
average effectiveness of about 0.6 in shear. 

* Grumman Aerospace Corporation. 
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INVESTIGATION OF THE EFFECTS OF PANEL IMPERFECTIONS 
ON A SIMPLIFIED MODEL 

Robert W. Herr and James C. Robinson 

Apparatus. - To gain some insight into the effects of panel imperfections on fuse- 
lage bending frequencies, the simplified model illustrated in figure 51 was fabricated. 

The six-bay model consists of a 0.51 mm (0.020 in.) aluminum panel specially selected 
for flatness, clamped between two frames of 1.02 by 19 mm (0.040 by 0.75 in.) alumi- 
num ai^le. The aspect ratio of each of the six bays (1.56) is representative of the aspect 
ratios of the 1/8-scale orbiter fuselage. Number 4 machine screws were used to clamp 
the panel between the frames and were spaced every 20 mm (0.79 in.) along the longerons 
and 19 mm (0.75 in.) along crosswise stiffeners. Panel imperfections were simulated by 
first loosening all of the clamping screws, then, with the frame lying on a flat surface, a 
distributed load was applied normal to the panel, one bay at a time, while the screws 
surroimding the bay were tightened. Upon removal of the static load, most of the static 
deflection remained. Adjacent bays were loaded in opposite directions. The nominal 
amount of deformation was dependent upon the magnitude of the distributed load. The 
procedure is approximate, of course, and the maximum simulated imperfection will vary 
somewhat from bay to bay; upon loosening of the machine screws, the panel always 
returned to its original flat state. 

For vibration tests, the panel model was oriented vertically with its upper end 
cantilevered. A 5 kg (0.0285 lb-sec2/in.) mass was clamped to the free end. To facili- 
tate the cantilever mounting and the attachment of the tip mass, the 0.51 mm (0.020 in.) 
panel extended 5.08 cm (2.0 in.) beyond the frame at each end. 

A small vibration exciter attached to the tip mass was used to excite the in-plane 
bending and the axial vibration modes. 

Analysis .- The effect of panel imperfections on the axial and bending stiffnesses 
of the structure shown in figure 51 was studied analytically usii^ the statics version of 
the Structural Network Analysis Program (SNAP /STATICS, ref. 38). One average lei^th 
bay of the beam measuring 38.1 cm deep by 2.44 cm long (15.0 by 9.625 in.) was modeled 
using triangular elements with both membrane and bending stiffnesses for the web and 
offset beam elements for the stiffeners. The model had a total of 160 grid points. A 
deformed shape of Wq[i - cos f(x) cos f(y^ was used for the panel imperfection. The 
resulting stiffnesses are for small displacements and low load levels inasmuch as the 
analysis does not consider the possibility of buckling. 

The effects of these stiffness reductions on natural modes and frequencies were 
calculated using the SPAR computer program (ref. 39) which is an improved version 
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of SNAP. A simple cantilevered-beam math model having 10 grid points was used to 
determine the lower natural frequencies and mode shapes. 

Results .- Panel imperfections ranging from 1.75 to 4 times the panel thickness 
caused the following computed reductions in the total stiffness (panel and stiffeners) of 
the simplified model: 


Panel imperfection 

Total stif 

fness reduction, 
percent 

Panel effectiveness, 
axial 

mm (in.) 

Wo/t 

Bending 

Shear 

Axial 

0.9 (0.035) 

1.75 

3 

6 

22 

0.63 

1.4 ( .055) 

2.75 

5 

12 

28 

.53 

2.0 ( .080) 

4.0 

8 

22 

32 

.48 


The panel imperfections investigated are seen to have only a minor effect on the model 
bending stiffness but a significant effect on the axial stiffness, with the effects on shear- 
ing stiffness falling in between. The effectiveness of an imperfect 0.51 mm (0.020 in.) 
panel (relative to a flat panel) in resistirg axial loads is given in the last column and was 
obtained by subtractirg the longeron stiffness from the computed axial stiffness of the 
complete model. The panel effectiveness is seen to range from 0.63 to 0.48 for panel 
imperfections between 1.75 and 4.0 times the panel thickness. 

The analytical and experimentally determined in-plane bending frequencies of the 
simplified panel model are as follows for the given values of initial panel imperfection: 


Panel imperfection 

Bending frequency, Hz 

Frequency reduction, 
percent 

Wo» 

mm (in.) 

^oA 

Experimental 

Analysis 

Experimental 

Analysis 

0 (0 ) 

0 

41.4 

44.7 

0 

0 

.9 ( .035) 

1.75 

40.3 

44.1 

2.6 

1.3 

1.4 ( .055) 

2.75 

39.2 

43.5 

5.3 

2.8 

2.0 ( .080) 

4.0 

38.6 

42.9 

6.8 

4.0 


The reduction in the in-plane bending frequency due to panel imperfection is seen to be 
quite small. This result is not surprising as the major out-of-plane panel imperfection 
lies near the neutral bending axis. 

For fuselage bending vibrations in the pitch plane, imperfections in the fuselage 
bottom panel would be expected to have a greater effect on the resonant frequencies 
since the entire panel area is in direct stress at a maximum distance from the neutral 
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axis. This effect is indicated in the following table in which the fundamental axial fre- 
quency of the simplified panel model is shown for various panel imperfections: 


Panel imperfection 

Axial frequency, Hz 

Frequency reduction, 
percent 

Wq, 

mm (in.) 

^o/t 

Experimental 

Analysis 

Experimental 

Analysis 

0 (0 ) 

0 

247 

279 

0 

0 

.9 ( .035) 

1.75 

196 

250 

21 

10 

1.4 ( .055) 

2.75 


240 


14 

2.0 ( .080) 

4.0 


232 


17 


The analytical results show a substantial reduction in the axial frequency for the range of 
panel imperfections investigated, 10 percent for imperfections 1.75 times the panel thick- 
ness to 17 percent for imperfections 4 times the panel thickness. Unfortunately, the 
experimental results are not definitive, due to coupling between the desired mode and the 
higher panel modes. In figure 52 the axial acceleration of the 5 kg (o.0285 Ib-sec^/in.) 
tip mass is plotted as a function of the excitation frequency for different amounts of panel 
imperfection. For the undeformed panel (fig. 52(a)) a relatively clean, well-defined reso- 
nance is observed at a frequency of 247 Hz. In figure 52(b) it is seen that for a nominal 
panel imperfection of 0.9 mm (0.035 in.) the frequency of the fundamental axial mode 
is 196 Hz, a reduction of 21 percent from the flat-panel frequency. For larger panel 
imperfections (figs. 52(c) and 52(d)), the out-of -plane vibration modes of the panels become 
so strongly coupled with the axial motion that it is impossible to isolate the desired axial 
vibration mode. 


The 21-percent reduction in the experimental axial frequency attributable to 0.9 mm 
(0.035 in.) imperfection corresponds to 37-percent reduction in the axial stiffness as com- 
pared with a computed 20-percent reduction in axial stiffness. If the longitudinal stiffeners 
are assumed 100-percent effective, the approximate effectiveness of the panel in carrying 
axial loads is given by 


Panel effectiveness = 



Al 


(E14) 


where the cross-sectional area of longerons, Al = 1.45 cm^ (o.224 in^), the cross- 
sectional area of panel, Ap 1.93 cm^ (o.3 in2), the fundamental axial frequency of 
flat panel model, fg = 247 Hz, and the fundamental axial frequency of imperfect panel 
model, f = 196 Hz. Substitution of these values into equation (E14) yields a panel effec- 
tiveness of only 0.34 compared with a computed value of 0.66. The cause of the large 
spread between the experimentally and analytically determined values of direct- stress 
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panel effectiveness is not immediately obvious, although there are reasons to believe that 
the true value probably lies between the two extremes. 

The 100-percent effectiveness of the longerons (or stiffeners) is a questionable 
assumption since the longerons do not connect directly to either the tip mass or the 
backstop. On the other hand, carry through of loads at the panel ends should be vei^ 
effective since axial extensions of the panel are clamped firmly to the tip mass and back- 
stop. This would suggest that most of the discrepancy between experimental and analyti- 
cal axial frequencies of the panel model with no panel imperfections (247 Hz and 279 Hz, 
respectively) is attributable to poor longeron carry -through structure. With the assump- 
tion of 100-percent effectiveness for the undeformed panel, the effectiveness of the 
longerons is given by 


Longeron effectiveness = 



(E15) 


where in this case f and fg are, respectively, the experimental and analytical axial 
frequencies of the panel model with no imperfections. Entering the appropriate values 
into equation (E15) results in a longeron effectiveness of 0.49. The effective-longeron 
area is thus reduced from 1.45 cm^ (o.224 in^) to 0.71 cm2 (o.ll in2)i Recalculation 
of the defer med-panel effectiveness using this reduced effective area of the longerons 
indicates an approximate panel effectiveness of 0.49. This value of panel effectiveness 
is in good agreement with results interpolated from the GAC nonlinear analysis (ref. 21) 
for a panel aspect ratio of 1.56 and an imperfection 1.75 times the panel thickness. The 
more sophisticated nonlinear analysis by GAC would be expected to predict the imperfect 
panel effectiveness more accurately than the linear SNAP analysis. 


CONCLUDING REMARKS 

Murray Bernstein* and John L. Sewall, Coordinators 

The estimate of direct stress capability based on 1/6 the web (shear panel) area 
lumped to the caps (longerons), which amounts to a 33 -percent-effective skin, yielded 
results that are too flexible as evidenced by model RSI. A more valid procedure by 
which the moments of inertia in the analytical model are reduced yielded model RS2. 
This method required the use of static test data. Model RS2 has a 46 -percent-effective 
skin fuselage and an 80-percent-effective skin wing in direct stress and correlates well 
with static test. 


* 


Grumman Aerospace Corporation. 
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Comparison of the rod and shear panel models RS4 and RS5 shows an increase in 
static deflection at the midfuselage with decreasing shear effectiveness. Vibration 
analyses of these models resulted in an increase in frequency with increasii^ shear 
effectiveness. Model RS4 exhibited frequencies closer to measured values but still 
somewhat higher. 

The membrane panel model RMl was determined to have an effectiveness of 
about 0.6 in both direct stress and shear. While model RMl clearly gives the best fre- 
quency and mode -shape correlation, the comparison with static test data is not so good 
as other static deflection models. 
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TABLE I.- PERTINENT SCALING RELATIONSHIPS FOR 


1/8-SCALE MODEL ^ 


Physical quantity 

Length, L, and displacements 

Poisson's ratio, v 

Mass density, p 

Modulus of elasticity, E . . 
Strain, e 

Area, A 

Area moment of inertia, I 


Magnitude 
Lm 1 

Lp ~8 

^m = ^p 
Pm - Pp 

^m _ , 
Ep 

^m = ^ 

Am ^ / 1\2 
Ap \8/ 



Mass moment of inertia, I' 
Stress, a 



Force, F ........ 

Longitudinal stiffness, EA 
Bending stiffness, El . . 
Torsional stiffness, GJ . 
Weight, W 



TeA)^ \8/ 

(EI)m /if 
(EI)p Uj 

(CrJ)m _ / 1\4 
(GJ)p - \ 8 ) 

Wp - is; 


Acceleration, a 

Natural frequency, o) 

^ Subscript m refers to model; subscript p refers to prototype. 


_ 8 

ap - 1 

‘^m 8 

o>p T 
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TABLE II.- MEASURED AND INITIAL ANALYTICAL MODE FREQUENCIES OF 

A 1/8-SCALE ORBITER MODEL 
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TABLE ni.- FIRST MODE NATURAL FREQUENCIES 
WITH STATIC PRELOAD 


Preload 

First mode 
frequency, 
Hz 

N 

lb 

0 

0 

43.7 

890 

200 

43.9 

1780 

400 

44.1 

2670 

600 

44.3 


TABLE IV.- ORBITER RUN STATISTICS 


(a) Static analysis 

Symmetric 


Field length, octal words 300 000 

CPU time for typical run, sec 2 800 

Number of runs 8 


Antisymmetric 
300 000 
2 750 
4 


(b) Dynamic analysis 


Phase I 
substructure 

Field length, 
octal words 

CPU time 
per run, 
sec 

Number of runs 

Symmetric 

Antisymmetric 

Fuselage 

250 000 

2680 

9 

1 

Wing 

160 000 

760 

6 

1 

Door 

160 000 

270 

2 

1 

Fin 

140 000 

160 

9 

2 

Payload 

160 000 

85 

1 

1 

Tape copy 

160 000 

60 

9 

1 

Phase H orbiter 

300 000 

3850 

9 

1 
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TABLE VI.- MEASURED AND BEST ANALYTICAL ORBITEH FREQUENCIES 


(a) Symmetric mode 


Mode 

Frequency, Hz 

Measured 

Analytical 

Model RMl 

Model RS6 

1 

43.6 

44.2 

44.2 

2 

51.6 

53.6 

54.4 

3 

57.9 

64.4 

63.0 

4 

1 

79.9 

84.2 

80.2 


(b) Antisymmetric mode 


Mode 

Frequency, Hz 

Measured 

Analytical 

Model RMl 

Model RS6 

1 

42.0 

44.5 

42.2 

2 

51.2 

54.3 

57.0 

3 

(a) 

67.4 

58.6 

4 

(a) 

77.8 

78.9 

5 

93.4 

88.7 

71.6 


^ Not found. 


TABLE VII.- ANTISYMMETRIC MODES VWTHOUT DOORS 


(a) Fuselage (NASTRAN, phase I) 


Mode 

Frequency, Hz 

Mode shape description 

Initial 

analysis 

(a) 

Model 

RMl 

Model 

RS6 

(a) 

1 

28.8 

24.4 

25.1 

First torsion, predominantly in the cabin 

2 

89.1 

75.4 

73.8 

First lateral bending 

3 

128.3 

111.7 

110.0 

Coupled lateral bending and torsion 


(b) Orbiter 


Mode 

Frequency, Hz 

Measured 

Analytical 

(b) 

1 

C27.0 

C23.7 

2 

(d) 

46.5 

3 

(d) 

57.5 

4 

(d) 

69.3 

5 

(d) 

79.6 

6 

(d) 

101.8 


^From reference 21, vols. 11, IHA, and IIIB. 
b Model RMl. 

c First torsion, predominantly in cabin. 

Not investigated. 







TABLE VIII.- DEFLECTION MEASUREMENT ACCURACY 


I 




Simultaneous loads applied (19 on-off cycles) 


Dial 

indicator 

Average 

deflection 

Standard 

deviation 

( 

K 

J -sigma deviation 

Ijm 

in. 

lira 

in. 

jum 

in. 

Percent 

1 

117.22 

4.615 X 10"3 

2.99 

0.118 X 10-3 

8.99 

0.354 X 10-3 

7.7 

2 

121.18 

4.771 

4.14 

.163 

12.40 

.488 

10.2 

3 

70.69 

2.783 

3.02 

.119 

9.04 

.356 

12.8 
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^ Not investigated. 

bUnits are jUm/Mm-N (in/in-lb). 




TABLE X.- STATIC DEFLECTIONS FOR SYMMETRIC LOADING 
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^Not investigated. 

bUnits are /im//im-N (in/in-lb). 





TABLE XIL- JOINT FLEXIBILITY EFFECTS ON 
ORBITER SYMMETRIC-MODE FREQUENCIES 


Mode 

Frequency, Hz 

Measured 

Analytical 

Initial analysis 

F2 

1 

43.6 

53.2 

51.9 

2 

51.6 

62.6 

63.4 

3 

57.9 

75.2 

72.55 

4 

79.9 

108.5 

91.8 




TABLE Xin.- PROPERTIES OF SIMPLIFIED FUSELAGE MODEL 


OF APPENDIX D 


tj 0.051 cm (0.02 in.) 

t2 0.0635 cm (0.025 in.) 

t0 0.051 cm (0.02 in.) 

b 0.3175 m (12.5 in.) 

0.419 m (16.5 in.) 

h]j 0.152 m (6 in.) 

Aj 0.645 cm^ (o.l in^) 

Ag 0.464 cm2 (o.072 in2) 

A 3 = A 4 0.564 cm2 (o.0875 in2) 

E 72.4 GN/m2 (l0.5 x I06 lb/in2) 

V 1/3 

i 2.46 m (97 in.) 


TABLE XIV.- FREQUENCIES OF SIMPLIFIED FUSELAGE MODEL 

OF APPENDIX D 




Frequency, Hz 

of - 


Mode 

Mathematical models 

NASTRAN model 

Bean 
solution 
(ref. 33) 

Piecewise 

Bending membrane 

Membrane, 


linear 

Coarse 

Refined 

coarse 

Symmetric, first 
vertical bending 

87.6 

77.8 

82.3 

78.5 

(a) 

Antisymmetric, first 
coupled transverse 
bending and torsion 

53.1 

51.8 

54.9 

51.1 

54.9 


^ Not investigated. 
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TABLE XV.- EFFECT ON PHASE II ORBITER SYMMETRIC-MODE 
FREQUENCIES OF REDUCED STIFFNESS 
IN FUSELAGE SIDEWALLS 



Frequency, Hz 

Mode 


Analytical 


Measured 

Initial analysis 

Minimal rod- 
shear panel 

1 

43.6 

53.2 

46.5 

2 

51.6 

62.6 

60.5 

3 

57.9 

75.2 

71.3 

4 

79.9 

108.5 

103.7 


3- Rods of minimal cross-sectional areas for panel stability. 




TABLE XVI. - EFFECTS ON FREQUENCY OF VARIATIONS 
IN PANEL EFFECTIVENESS 


(a) Symmetric mode 



Frequency, Hz 

Mode 

Measured 

Initial 

analysis 

RS3 

RS4 

RS5 

RS6 

(a) 

RMl 

1 

43.6 

53.2 

53.85 

45.55 

46.8 

44.2 

44.2 

2 

51.6 

62.6 

59.5 

54.3 

55.1 

54.4 

53.6 

3 

57.9 

75.2 

74.4 

63.45 

65.8 

63.0 

64.4 

4 

79.9 

108.5 

108.7 

86.3 

88.3 

80,2 

84.2 


(b) Antisymmetric mode 


Mode 

Frequency, Hz 

Measured 

Initial 

analysis 

(a) 

RMl 

RM2 

RS6 

(a) 

1 

42.0 

52.9 

44.5 

44.4 

42.2 

2 

51.2 

72.6 

54.3 

54.2 

57.0 

3 

(b) 

85.1 

67.4 

67.2 

58.6 

4 

(b) 

92.0 

77.8 

76.5 

78.9 

5 

93.4 

101.5 

88.7 

86.4 

71.6 


^From reference 21. 
b Not found. 
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TABLE XVn.- EFFECTIVE CROSS-SECTIONAL AREAS FOR PANELS IN PAYLOAD BAY 


O 

C 
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^•Area of sidewall, Dts; area of lower skin panel, Bt^ 






85 


Figure 1.- Photographs of 1/8 -scale dynamic model of space -shuttle configuration 
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2.- Photograph of 1/8-scale orbiter model with cargo 
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Figure 3.- Photograph of l/8-scale orbiter model with cargo-bay door installed. 
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(139.5 in.) 

Figure 4. - General arrangement drawii^ of 1/8-scale orbiter model. Dimensions are normalized 

to fuselage overall length 3.543 m (139.5 in.). 
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Figure 5.- Construction of l/8-scale orbiter fuselage. 
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Figure 6.- Strain-gage installation on skin panels and loi^erons. 






L-73-26411 

Figure 7.- Partially constructed wings of 1/8-scale orbiter. 
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F^ure 8.- Setup for vibration tests in horizontal attitude. 
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Figure 10.- Schematic drawing of fixed transducer locations for vibration tests. 








Figure 11.- Schematic diagram of instrumentation. 
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97 


Measured and analytical mode shapes tor the orblter second symmetric mode. 
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Figure 14.- Measured and analytical mode shapes for the orbiter third symmetric mode. 
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Figure 15.- Measured and analytical mode shapes for the orbiter fourth symmetric mode. 




strain. 








strain 



Fuselage 

loi^erons 

Gages (LH), sta. 116.0 

O RH side 
□ LH side 

Door 

J8 longerons 


Load. N 


tllu 


Load, lb 

Figure 17.- Static strain-gage measurements obtained on cargo-bay door 

and fuselage longerons.- 






Ballast 



Normalized amplitude 

Figure 18.- Effect of fin ballast on experimental first- mode frequency and motion. 
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(a) Main substructures. 

Figure 19.- 1/8-scale orbiter NASTRAN model (NASTRAN element identifications in parentheses). 
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support brackets assumed rigid). 
Figure 19.- Concluded. 






Frames in engine support structure 

Figure 20.- Typical NASTRAN rod-shear panel models for frames 

for 1/8 -scale fuselage. 
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Total central processing unit time for full calculation 7915 sec. 

Figure 21.- Orbiter substructurii^ sequence and computation times. 
(Phase in not executed in this investigation.) 
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(a) Overall view into fuselage interior. 

Figure 22.- Views of wing carry -through structure with wing removed. 
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(b) Wii^ -fuselage interface at fuselage station 144.75. 
Figure 22.- Concluded. 






(a) Forward and center spar connections. 


Figure 23.- Fin-fuselage interface. 
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L-74-3989.1 

(b) Detail of aft fin-fuselage interface connection. 

Figure 23.- Concluded. 
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Figure 25.- Model of cabin-to -cargo -bay joint flexibility. 
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Figure 26.- Measured and analytical mode shapes of model RMl for the orbiter first symmetric mode. 
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Figure 27.- Measured and analytical mode shapes of model RMl for the orbiter second symmetric mode. 
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Figure 29.- Measured and analytical mode shapes of model RMl for the orbiter fourth symmetric mode. 
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Figure 30.- Measured and analytical mode shapes of model RMl for the orbiter first antisymmetric mode. 
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(b) Vertical motion (z-direction). 
Figure 32.- Concluded. 
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Figure 33.- Measured and analytical mode shapes of model RMl for the orbiter first 
torsional mode with the cargo door removed. 
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Figure 34.- Support of 1/8-scale orbiter model for static-load tests. 
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(a) Forward fuselage connection. 

Figure 35.- Orbiter-model interstage fittings. 
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(b) Aft fuselage connection. 
Figure 35.- Continued. 
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Figure 35.- Concluded. 
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Figxire 36.- Schematic of orbiter support point constraints for static tests 



127 


(a) Symmetric. 

Figure 37.- Schematic of load application points and loading directions. (See tables IX, X, and XI.) 
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Figure 37.- Concluded. 
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(a) Forward fuselage. 

Figure 38.- Typical loading apparatus and dial -indicator setup for static -load tests, 
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Figure 38.- Concluded. 
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Figure 39.- Schematic of basic dial -indicator locations for deflection measurement. 














Load, lb 














Midspan 
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Figure 42.- Local cross-sectional coordinate system of appendix D. 





Kfeel Floor Sidewrdi Keel Floor Sidewall 

(a) Coarse model. (b) Refined model. 

Figure 43.- NASTRAN models of simplified fuselage of appendix D. Dimensions are in m (in.). 
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Figure 44.- First symmetric lateral bending mode for NASTRAN (piecewise linear) simplified fuselage 

of appendix D; f = 78.5 (77.8) Hz. 
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Figure 46.- NASTRAN simplified 1/4 -model fuselage of appendix D with 
high-aspect-ratio QDMEM2 elements. Dimensions are in m (in.). 
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Aft fuselage and 



Fuselage station 

Figure 47.- Comparison of static measured and analytical deflections due to vertical loading 

at fuselage nose. 
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Wing vertical deflection. 



0 10 20 30 40 50 60 ii 

Wir^ buttock line at station 162 

Figure 49. - Comparisons of static measured and analytical deflections 

for wir^ tip loadii^. 



Total load, 


Lengthening 


Shortening 




Compression 


^ 1 ' 




■'/y/ 


0 .1 .2 .3 


0 20 40 60 80 

End load 


200 160 120 80 40 0 40 80 120 160 200 

End deformation e , pm/m 

(a) Fuselage bottom surface; typical aspect ratio a/b = 1.32; 
t = 0.635 mm (0.025 in.). 

Figure 50.- Analytical panel effectiveness estimates from reference 21, vol. II; 

E = 68.95 GN/m2; v = 0.3. 




Lengthening ^ Shortening 



End deformation e , jum/m 

(b) Fuselage sidewall; typical aspect ratio a/b = 1.79; t = 0.635 mm (0.025 in.). 

Figure 50.- Continued. 













Figure 51.- Simplified panel model for panel effectiveness studies 







Figure 52.- Longitudinal response of tip mass as function of excitation frequency for 
various out-of-plane static panel deflections. 






Frequency, Hz 

(d) Panel deflection of 1.8 mm (0.071 in.). 
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